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ABSTRACT: Offshore gravity structure foundations, including their analysis and 
design, are investigated. The results experimental study the interaction 
such structures with their soil foundations when subjected wave forces are reported. 
The experimental tests were performed centrifuge 100 centrifugal acceleration 
with 10-cm diam model structure bearing normally consolidated desposit 
Santa Barbara silt. Therefore, the experimental test results apply 10-m diam 
prototype structure. The model structure was subjected both monotonic and cyclic, 
well inclined and eccentric loading conditions; and its resulting movements 
(sinking, sliding and tilting) were recorded. 
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ABSTRACT: Experimental and analytical studies are described the offshore gravity 
structure foundations with their soil foundations when subjected wave forces. The 
experimental tests were performed centrifuge 100g centrifugal acceleration with 
10-cm diam model structure bearing normally consolidated deposit Santa 


Barbara silt, and therfore apply 10-m diam prototype structure. The analytical 
predicitions, for the particular prototype situation considered this study, were 
obtained using finite element program which incorporates three-dimensional 
nonlinear soil model capable describing the behavior the Santa Barbara silt under 
both monotonic and cyclic loading conditions. The results both two-dimensional 
(plane strain) and three-dimensional calculations are reported, and the are compared 
with the experimental test results obtained the centrifuge. 
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ABSTRACT: The generation and propagation vibrations due construction are 
discussed. The relative intensity the vibration produced many construction 
sources identified. brief summary the development suitable vibration- 
measuring instruments related. Differences current damage criteria are revealed. 
Future research needed the effects both transient and steady-state vibrations 
structures, uncured concrete, soils and underground utilities; additional research 
needed fatigue building materials. Human response evaluation transient 
vibrations also requires further research. 


REFERENCE: Wiss, John F., “Construction Vibrations: State-of-the-Art,” Journal 
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KEY WORDS: Dynamics; Finite element method; Pipelines; Seismic 
design; Soil-structure interaction; Tunnels; Wave propagation; Waves 


ABSTRACT: The response buried long structures waves traveling the 
longitudinal direction the structures was studied. The structures and the surrounding 
soils were represented three-dimensional finite elements. The three-dimensional 
models, however, were reduced two-dimensional models using special type 
wave-transmitting element which faithfully accounts for the effects wave 
propagation. The results obtained for linear system indicate that, for large structures, 
soil-structure interaction effects the computed stresses and strains are indeed 
significant; therefore, will overly conservative use free field strains the 
design. Also, the assumption vertically propagating waves not adequate 
predicting response; therefore, appropriate considerations should given traveling 
waves. 


REFERENCE: Nwang, Richard N., and Lysmer, John, “Response Buried Structures 
Traveling Journal the Geotechnical Engineering Division, ASCE, Vol. 
107, No. GT2, Proc. Paper 16052, February, 1981, pp. 183-200 


16063 AT-REST PRESSURE PEAT SOILS 


KEY WORDS: pressure; Coefficient earth pressure; 
Compressibility (soils); Fibers; Lateral pressure; Peat; Soil classification; Soil 
compression tests 


ABSTRACT: The results are presented experimental program directed toward 
evaluating (the coefficient lateral earth pressure under conditions lateral 
formation) and its dependence stress history and other parameters peat soils. 
Three peats with different fiber contents were consolidated custom-built Ko-test 
tube. During loading and rebound, the axial deformations, pore pressure, and both axial 


and radial stresses were measured. sand and clay sample compared well with 
previously reported values and empirical predictions. The values for the peat 
sample were found remain essentially constant during the loading phase, but their 
values (which increased gradually during unloading) never exceeded unity. Several 
indirect approaches were tried compute the the peats from effective shear 
strength parameters. Peat type appeared exert significant influence the 
measured values during loading. 


REFERENCE: Edil, Tuncer B., and Dhowian, Abdulmohsin W., “At-Rest Lateral 
Pressure Peat Soils,” Journal the Geotechnical Engineering Division, ASCE, Vol. 
107, No. GT2, Proc. Paper 16063, February, 1981, pp. 201-217 
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INTRODUCTION 


The fast-growing need for new oil and gas supplies has led the development 
numerous new offshore fields. Until the last decade most offshore petroleum 
industry activity was confined areas shallow water depths. few years 
ago active offshore drilling and production began the North Sea and now 
being contemplated both the North Atlantic and Arctic waters. these 
sites the waves are larger and the water depths are greater than areas previously 
exploited. New types offshore platforms have thus had designed 
order withstand such difficult site conditions. Some these structures are 
unprecedented dimensions, and many the problems encountered their 
design have never been tackled before. particular, the design their 
foundations poses new and very complicated geotechnical problem, explained 
later. 

Prior North Sea oilfield development, the majority offshore structures 
has been tubular steel design, pinned the sea floor piles, e.g., the 
Gulf Mexico. Initially, platform structures developed for the North Sea were 
also this type, but substantially stronger construction than those built 
the Gulf Mexico. The soil conditions found some other parts the 
offshore world are different general from those the Gulf Mexico, and 
often from foundation point view are excellent that, encountered 
land, would not have required piles for the foundations even heavy 
structures. Offshore, the piles are needed because the very large lateral 
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and overturning loads generated waves. The stability such platforms depends 
the presence the piles. However, extremely difficult and time consuming 
drive piles deep and rough waters. 

Consequently, there potentially hazardous period regions where bad 
weather prevalent, after the placement the platform, and before sufficient 
piles are driven ensure the safety the structure. This period may months 
long. There were, therefore, good reasons consider alternative design solutions 
for major offshore structures the North Sea—designs for which the structure 
would stable within days after installation. Studies demonstrated that 
many locations should possible place gravity-type concrete structures 
directly the sea floor without the use piles. number such concrete 
structures have now been designed and built. 1978, concrete gravity 
platforms had been installed the North Sea water depths ranging from 
100 m-150 The structures are built shore, floated the site, sunk and 
seated the sea floor ballasting. They are designed resist wave forces 


FIG. 1.—Gravity Platform Load Conditions and Definitions Terms 


the combined action their weight and large foundation dimensions (typically 
the order several hundred thousand tons and 10,000 respectively). 

terms size these gravity structures are one order magnitude larger 
than platforms previously built, and this also applies the environmental forces 
which they are subjected. Fig. example, gives the order magnitude 
result its interaction with wave whose height and period sec. 
When such wave passes the structure, the forces first act one direction 
and few seconds later the opposite direction with almost the same magnitudes. 
Under such cyclic loading conditions, soils exhibit very significant dissipative 
and softening characteristics (1,3). The effect cyclic structural forces 
soil foundations thus type fatigue problem, and great concern 
determine what extent the behavior soils under cyclic loading may 
influence the safety these structures. 

Field data the performance full scale structures under severe loading 
conditions are urgently needed, but until now have remained proprietary. Further 


000 
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data are unlikely obtained soon due the scale the structures involved, 
the cost testing, and the low probability having particular instrumented 
structure subjected design loads. Consequently, some form model study 
seems inevitable enable alternative design analysis methods checked. 

program whereby this possibility could implemented was suggested 
Chevron Oil Company. would utilize centrifugal model tests gravity 
platform structure placed normally consolidated deposit soil carried 
out California Institute Technology. The centrifuge model tests would 
include both monotonic and cyclic loading the model structure. Centrifuge 
tests were proposed because them, model scaling relations are adequately 
satisfied considered later. The soil involved was clayey silt, obtained from 
the ocean floor off the southern California coast. Centrifuge tests the behavior 
gravity platform placed homogenous clay deposit have been reported 
before Ref. 

This proposal was accepted Chevron Oil Company, and work was initiated. 
This report describes the apparatus, the tests and test results. Their analysis 


the scaling relations required characterize tests geotechnical models 
are established, found that, because the general dependence the 
mechanical properties soil the ambient stress conditions, and the importance 
gravity-induced stresses, scaling can only satisfied under special conditions. 


the special case soil, inconvenient impossible construct model 
material, and real soil usually employed model tests. that case, the 
scaling conditions require that the soil model subjected higher gravitational 
acceleration than the prototype. The ratio the accelerations model and 
prototype structures inversely proportional the ratio their linear dimen- 
sions; that say, 1/100-scale soil model must subjected 100 times 
the earth’s gravitational acceleration, and on. obtain the necessary 
accelerations, centrifuge required. Considering the size and mass soil 
model experiments, even 1/100 scale, the centrifuge has quite large. 

the ratio linear prototype dimensions those the centrifuge model 
then the ratio area and volumes The scaling relations indicate 
that forces the prototype are times those the model, that stresses 
(force per unit area) are unchanged. Deformation the prototype times 
larger than the model, but strains (deformation per unit length) are the same. 
Thus, the presence the same material both prototype and model results 
identical stresses and strains homologous points. Where dynamic problems 
are involved, turns out that time the prototype times the time 
the model. Consequently, model frequencies are higher the factor but 
velocities are unchanged. Energy the prototype times the energy 
the model but energy density (energy per unit volume) the same. Table 
gives more complete list the relations between prototype and model 
(centrifuge) parameters. 

The attractiveness the centrifugal method that the stresses the model 
are identical those the prototype that avoids problems associated 
with testing, earth gravity, small soil models involving material with strongly 
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nonlinear and stress-dependent behavior. The disadvantages are associated with 
performing the tests models which are rotating rates 100 
rpm centrifuge. Power and signals have passed and out through 
electrical and hydraulic sliprings. till now, electrical noise has presented 
problem recording the results variety experiments, especially where 
strain gauges have been employed. The noise comes from ambient sources, 
the electric motor driving the centrifuge, from motors used control the tests, 
and centrifuge vibration. 

number centrifuges have been built and used for soil testing. There 
are three the United Kingdom, two Cambridge and one Manchester, 
with radii and acceleration capabilities 200 the Soviet 
Union, recent paper (4) refers the employment centrifuges 
for soil testing purposes. far, only one two small centrifuges have been 


TABLE 1.—Scaling Relations 


Full scale 
(prototype) 
(2) 


Centrifugal 
model, gravities 


(3) 


Quantity 
(1) 

Linear dimension, displacement 

Area 

Volume 

Stress 

Strain 

Force 

Mass 

Acceleration 

Energy 

Density 

Energy density 

Velocity 

Time: dynamic terms 

Time: diffusion cases 

Time: viscous flow cases 

Frequency dynamic problems 


used for such tests the United States, although the technique was apparently 
originated here (2); the reasons for this limited usage have not been determined. 
Centrifuges are currently used for geotechnical research Sweden, Denmark, 
France, and Japan, apart from the Soviet Union and the United Kingdom. 
compilation references centrifugal testing, worldwide, extends more 
than 150 papers and number books. 


Testeo 


The soil tested this investigation sandy clayey silt Barbara 
obtained soil investigations off the California coast the vicinity 
Santa Barbara and supplied Dames and Moore, San Francisco, through the 
courtesy Chevron Oil Research Company, Habra, Calif. The scale tests 
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were carried out the soil rectangular container 14.5 in. 11.5 in. 
in. mounted the centrifuge platform. The Santa Barbara silt initially 
obtained, was supplied the form dried out chunks soil, left over from 
laboratory tests performed Dames and Moore. Later, extra quantities 
this material were received moist condition plastic sample tubes, from 
which the soil was moved jacking with piston extractor. was not convenient 
correct either case, place the soil the test containers the form 
which was received, since was desirable install the platform 
normally consolidated soil deposit simulate ocean floor test site. Conse- 
quently, both dry and wet as-received samples Santa Barbara silt were mixed 


depth 
soil model (cm) 


FIG. Shear Strength Profiles Series Tests 


with tapwater mixing machine until the resulting material was smooth and 
lump-free water content slightly higher than the liquid limit the soil 
(w, 30%). 

The soil slurry was then poured into the test container, taking care avoid 
the inclusion air bubbles. prepare the test specimen the container was 
centrifuged 100 the test acceleration, until had consolidated. Drainage 
was facilitated the presence sand layer the bottom the container, 
overlaid blotting paper prevent mixing the silt and sand, and 
layer porous plastic placed inside the wall the test container. 

The question the similarity the centrifuge test specimen soil that 
the larger prototype section, and particular, the identification the stress 


7 
——T5 
10. 
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paths model and prototype complicated. Absolutely similar stress histories 
and stress paths cannot obtained during testing, but reasonable facsimile 
can achieved the centrifuge (8). 

The soil the site slightly overconsolidated with field shear strength 
the line about 2.5 The strength increases linearly with depth 
The field strengths were measured performing tests the samples 
they were obtained the deck the drilling vessel. such, they are 
not very precise and depend largely the local nature the soil exposed 
the end each sample tube for testing. 

order effect close simulation the real material properties the 
field, the soil the container was slightly overconsolidated covering with 
thick layer sand during the consolidation phase. The sand was removed 
before performing the model tests. 

The duration centrifuging was arrived by: (1) Consideration the 
consolidation properties the soil obtained from consolidation tests; and (2) 
performing miniature vane shear and water content tests samples from 
varying depths the test container. When the soil fully consolidated, the 
shear strength profile should increase essentially linearly with depth, and the 
water content should decrease with depth. Both dry and total unit weight should 
increase with depth. 

One set soil test results from the containers which the model tests 
were performed shown Fig. The silt occupied the container depth 
12.50 cm, representing depth 12.50 the prototype. The soil was 


probed about from the sides the container. apparent that 


BASKET ASSEMBLY 


SOIL BASKETS 
HEIGHT, 10° 


wom 


FIG. 3.—Perspective View Centrifuge Test Container 
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most cases adequate consolidation has been achieved although the soil strength 
appears slightly lower than the field. This appears related 
the changing soil-water chemistry caused continually diluting the original 
salt content the soil the addition tapwater. 


AND INSTRUMENTATION 


Centrifuge.—The centrifuge used Model A1030 Genisco 
which consists 80-in. diam aluminum alloy arm which rotates the horizontal 


FIG. 4.—Centrifuge Model Platform and Loading Devices 


plane and rated 10,000 g-lb payload capacity. each end the arm 
located in. in. magnesium mounting basket capable carrying 
payload 100 60-lb payload 175 The acceleration range 
the approx 40-in. radius the baskets from variable speed 
transmission based the Rouverol ball-galaxy principle allows continuous 
variation the acceleration throughout the range. The transmission driven 
1,800 rpm 220 three-phase double-ended electric drive motor. 
For accurate determination the rotational speed, 600 tooth gear wheel 
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located the main drive shaft and this produces 600 electrical pulses per 
revolution magnetic pickoff. These pulses are read electronic counter 
which converts them LED display RPM accurate 0.1 prm. The 
drift and wow the system any given setting 0.05%. The acceleration 


FIG. 6.—Definition Test Parameters 


arm housed extruded aluminum enclosure, with all controls and 
instrumentation, the interests safety, located remotely. 

Electrical power and signals and from the rotating arm basket are conducted 
through sliprings, various capacities the amp-30 amp range. Externally 
generated hydraulic air pressure transmitted through either two four 


FIG. 5.—Centrifuge Model Platform 
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lines means rotary unions (hydraulic sliprings). Operations the centrifuge 
can observed means television camera mounted the rotating 
arm close the axis; its signal conveyed either through the aforementioned 
sliprings through coaxial cable and related, separate sliprings, monitor 
television the instrumentation house. 

The centrifuge layout shown Fig. with more detailed view the 
arm the test configuration this report Fig. which will described 
later. 

Specialized Test Apparatuses.—Fig. shows one the two apparatuses which 
were constructed apply inclined and eccentric force the model platform. 
Loading was achieved directly pushing the model platform with small 
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FIG. T1: Test Results 


hydraulic piston attached the aluminum loading frame shown Fig. This 
loading frame bolted the aluminum container. 

The model platform was machined out aluminum 10-cm diam shown 
Fig. Its weight was adjusted that 100 corresponded bearing 
pressure 4.3 The platform was mounted with skirts (0.4 thick) 
whose depth could varied cm-0.6 cm) adjusting the position the 
movable base the structure. The bearing surface consisted aluminum 
plate which sand had been glued that would rough simulate 
field conditions. 

The applied force and the resulting movements the model platform were 
recorded load cell and displacement transducers. The cell consists 
ring stainless steel ring’’), 1.0-in. 0.5-in. height 0.025-in. wall 
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thickness, which were bonded four strain gauges (BLH FDE-06S-35ET, G.F. 
3.10 0.5 ohms, Serial No. 1-A-KH, Lot No. A-D-135) using epoxy 
adhesive (BLH EPY-150). Two the gauges were located the outside 
the ring diametrically opposite one another, while the other two gauges were 
placed the same location, but the inside. This arrangement provided 
reasonable output and good temperature compensation when wired into complete 
bridge. Displacements were recorded through nylon strings attached the 
platform and connected displacement transducers made small cantilever 
strip metal which strain gauges had been glued. They formed the legs 
Wheatstone bridge. The displacement transducers were preloaded slightly 
avoid any slack the strings. 
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FIG. 8.—Test T5: Test Results for Vertical Load Versus Vertical Displacement 


Data Acquisition Systems.—The recording instruments used were: (1) Hon- 
eywell Model 1858 CRT visicorder, which allows inertialess recording from 
5,000 Hz; and (2) Hewlett-Packard Model 7045A X-Y recorder. The 
X-Y recorder was used primarily check that the tests were proceeding 
expected; accordingly only the load and one displacement component were 
recorded it. The load and all displacements were also recorded the visicorder. 

All pretest calibrations were carried out using the entire electronic circuitry, 
i.e., the calibration signals were routed through those terminals, amplifier 
channels, and sliprings which they would use during the actual testing. calibrate 
the load cell, brass tester weights were suspended from the cell and the output 
recorded both the X-Y recorder and visicorder. The calibration the 
displacement transducers was accomplished with the aid Federal dial gauge 
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accurate 0.001 in. The output was again recorded both the X-Y recorder 
and the visicorder. 

Test Procedures.—The Santa Barbara soil was thoroughly mixed with water, 
previously described, and placed the prepared soil container. The preparation 
consisted placing saturated sand the bottom the container 
and covering both the sand and the vertical walls the cylinder with porous, 
permeable plastic material. This facilitated the drainage water 100 
maximize consolidation. Once the soil was introduced the centrifuge was brought 
100 acceleration and allowed operate this setting for several hours 
consolidate the soil. The centrifuge was then stopped, and for the initial 
tests vertical marker columns wax were installed the middle plane 
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FIG. 9.—Test T5: Test Results for Moment Versus Tilt 


the bucket pouring hot wax predrilled holes about apart. Later, 
modeling clay cylinders were used place the wax. The soil deposit was 
thereafter covered with 2-1/2-cm layer sand top blotting paper, 
and the centrifuge was again brought 100 acceleration for several hours. 
The consolidation duration was established from consolidation tests the clay 
and trial runs the centrifuge with shear strength measurements. After consoli- 
dation, the centrifuge was stopped, the sand layer removed, the aluminum loading 
frame bolted the top the container, and the model platform installed 
top the silt deposit. Electrical connections were established and the centrifuge 
returned 100 acceleration. The system was allowed stabilize this 
acceleration for few minutes, during which time the recording instruments 
were allowed warm up, and the calibrated excitation voltages were set (typically: 
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load cell 6.00 d-c; displacement transducers 1.57 d-c). The centrifuge 
was then stopped again, and the loading device was set the desired inclination 
and eccentricity the applied force. The displacement transducers were 
connected the platform and the centrifuge brought again 100 acceleration 
ready for testing. The tests are described the next section. 

Data Reduction.—The strip chart from the visicorder was digitized onto punched 
cards (Benson-Lehner 099 data reducer 282E Telecordex IBM cardpunch) 
which were used with computer program yield various plots load versus 
displacement the model platform. After test had been completed, the 
centrifuge was stopped, the loading frame dismounted, and the soil deposit 
cut half through the plane the wax modeling clay columns investigate 
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FIG. 10.—Test T6: Test Results for Vertical Load Versus Vertical Displacement 


the resulting movements the soil foundation. Miniature Tor-vane, water content, 
and dry density tests were conducted various depths the soil deposit. 


total tests (6) were performed various inclinations and eccentricities 
the applied external (wave, ice) force Fig. shows the notation and 
defines the parameters used the interpretation these tests. The load 
conveniently identified means its horizontal and vertical components, 
and respectively, giving rise the moment Under the action 
the load the platform sinks, slides and tilts, and and [in Fig. 
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FIG. 12.—Test T6: Test Results for Moment Versus Tilt 
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FIG. 13.—Test T6: Test Result for Distorted Configuration Soil Deposit Cross 
Section after Failure 
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FIG. 14.—Test T8: Test Results for Vertical Load Versus Vertical Displacement 
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FIG. 15.—Test T8: Test Results for Horizontal Load Versus Horizontal Displacement 
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are used denote its resulting vertical, horizontal and tilting displacements. 
The test results are grouped Ref. and presented the form load versus 
displacement curves for the various deformation modes. Also included are the 
soil profiles corresponding each test. illustration, Figs. 7-13 show 
and 10°; —11%), respectively. tests through the platform 
was first subjected few cycles small magnitude loading before the force 
was monotonically increased until failure the platform occurred, and, 
illustration, Fig. shows the distorted configuration the soil foundation 
test failure and clearly indicates the existence failure surface. 
test 10°; +11%), the model platform was first subjected 
cycles constant magnitude loading then six further cycles loading 
with larger magnitude before the load was monotonically increased failure 
(Figs. 14-16). From Figs. apparent that upon cycling the force 
the structure accumulates permanent vertical (Fig. 14), horizontal (Fig. 15), 
and tilting (Fig. 16) movements and that the rate accumulation decreases 
the number cycles loading increases. After six cycles loading the 
rate accumulation irreversible displacements very small and eventually 
would vanish further cycles loading with the same magnitude were applied 
the platform. This the phenomenon referred 


This investigation concerns offshore gravity structure foundations, their analy- 
sis and design. The results experimental study the interaction such 
with their soil foundation when subjected wave (or ice) forces 
are reported. The experimental tests were performed centrifuge 100 
centrifugal acceleration with 10-cm diam model structure bearing normally- 
consolidated deposit Santa Barbara silt. The experimental test results therefore 
apply 10-m diam prototype structure. The model structure was subjected 
both static and cyclic, inclined and eccentric loading conditions, and its resulting 
movements (sinking, sliding and tilting) were recorded. proved possible 
carry out the investigation the small centrifuge California Institute 
Technology, adequately simulating the response prototype 10-m diam platform 
soil prototype-scale properties. Although this not the size the prototype 
production platforms the North Sea, sufficiently large not considered 
model test laboratory scale, and the results are useful their own right. 
Since the tests can set and run safely, relatively quickly, and reliably 
the small centrifuge, its use further experiments aid assessing 
the reliability alternative design analysis techniques encouraging. 
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OFFSHORE GRAVITY STRUCTURES: ANALYSIS 


ASCE, and Ronald ASCE 


INTRODUCTION 


terms size, offshore gravity structures are one order magnitude larger 
than platforms previously built, and this also applies the environmental forces 
which they are subjected. When wave passes such structure, the forces 
first act one direction and few seconds later the opposite direction 
with almost the same magnitudes. therefore not sufficient check the 
stability the structure for the maximum forces introduced say, 100-yr 
wave acting static load, but must also considered that, before the 
maximum wave occurs, the structure has already been subjected great 
number (several thousand) lesser waves causing cyclic loading the soil 
foundation. Under such cyclic loading conditions, soils exhibit very significant 
dissipative and softening characteristics (1,9). The effect cyclic structural 
forces soil foundations thus type fatigue problem, and great 
concern determine what extent the behavior soils under cyclic loading 
may influence the safety these structures. 

order that the design the foundation rational the interaction between 
the structure and its soil foundation under cyclic, inclined, and eccentric loading 
conditions must analyzed. The complexity the problem dictates that the 
method analysis numerical nature. 

make accurate predictions, realistic analytical model soil stress-strain- 
strength behavior essential. Indeed, must emphasized that the solution 
provided even the most complex finite element program will never better 
more realistic than the constitutive equations adopted the program 
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model soil behavior. particular, the cyclic wave forces subject the foundation 
materials three-dimensional stress-strain histories involving cycles loading, 
unloading, and subsequent reloading. Consequently, for the analysis the 
foundation-structure interaction realistic, the dependence the constitutive 
equations their stress-strain history must taken into account. However, 
most the stress-strain strength soil models that have been used far 
geotechnical engineering cannot cope adequately with such complex loading 
histories. Various simplified methods analysis have thus been proposed and 
are currently use. 

The first category methods based failure concepts. common 
practice soil engineering design structure foundations with adequate 
factor safety against which general sense must interpreted 
the structure’s loss equilibrium with uncontrolled displacements. For 
gravity structure, three modes failure can identified, namely, bearing, 
sliding, and overturning. The corresponding failure loads are computed 
assuming that the actual three-dimensional field situation can idealized 
two-dimensional (plane strain) representation, that the soil foundation behavior 


a 
10 


FIG. 1.—Santa Barbara Silt: Compression and Axial Extension Tests Results 


can taken isotropic and, rigid-perfectly plastic, and using limit 
equilibrium techniques (see e.g., Refs. 11, 12, and 26). The results 
the failure analysis are usually presented the form interaction diagram 
(26). The soil strength used interpreting such diagram for design 
purposes strength account for the effects cyclic loading. 
This method analysis evidently very crude and its results are rather uncertain 
since requires many simplifications the actual field conditions. 
therefore not clear what extent analytical predictions based such drastic 
simplifying assumptions are related actual field performances. Furthermore, 
the method does not provide information about deformations prior failure. 
attempt provide this, more complicated analytical methods have 
been proposed based finite element representations the soil foundation 
underneath the structure. the studies, soil modeled either isotropic, 
linear, nonlinear (3,5) elastic material, isotropic, elastic-perfectly plastic 
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(25) elastic-plastic hardening (27) material. order account for the effects 
cyclic loading the soil foundation, soil specimens the laboratory are 
subjected the cyclic shear stresses predicted the finite element analysis, 
and the results these tests are used predict the resulting deformations 
the platform foundation. This method analyses therefore essence 
procedure (8). Although the stress path method may have proven 
useful soil engineering, its limitations should recognized, namely that the 
estimation the stress (or strain) path followed any soil element the 
gound strongly dependent the particular soil model used for these predictions. 
This stress path then employed test the soil determine its strain 
response, which usually widely different from that used the stress calculation. 
therefore evident that the meaning such procedure best limited. 
Recently more rational analysis procedure has been proposed Refs. 
and 20. The method based finite element analysis the soil-structure 


TABLE 1.—Santa Barbara Silt: Yield Surface Parameters, 133.333 


(m) , 


(4) 


4 


interaction problem using soil model which describes the nonlinear, 
anisotropic, elastoplastic, path-dependent stress strain-strength properties soil 
when subjected complicated, and particular cyclic, three-dimensional 
loading paths. This formation allows the solution the complicated wave-struc- 
ture-soil interaction obtained directly because incorporates the material 
stress-strain-strength properties, including their progressive softening with strain. 
However, although this formulation has been carefully checked against laboratory 
experimental test results, its accuracy predicting field performances has not 
yet been tested. This common problem geotechnical engineering. exists 
because the paucity field test information. Field data the performance 
full-scale structures under severe loading conditions are urgently needed, 
but until now have remained proprietary. Further data are unlikely obtained 
soon because the scale the structures involved, the cost testing, and 
the low probability having particular instrumented structure subjected 
design loads. Consequently, some form model study seems inevitable 
enable alternative design analysis methods checked. 

program whereby this possibility could implemented was suggested 
Chevron Oil Company. would utilize: (1) Some the results laboratory 


(1) (2) (3) (5) (6) 
0.167 0.0 0.000 266.667 
0.170 0.0 0.090 69.565 
0.152 0.0 0.205 37.037 
0.140 0.0 0.270 21.739 
0.130 0.0 0.330 13.605 
0.120 0.0 0.360 9.697 
0.123 0.0 0.393 6.617 
0.115 0.0 0.415 4.737 
0.117 0.0 0.465 2.614 
0.123 0.0 0.502 1.821 
0.134 0.0 0.550 0.0 
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tests particular soil performed University California Los Angeles 
(UCLA) (9); (2) centrifugal model tests model structure placed normally 
consolidated deposit the same soil carried out California Institute 
Technology (17,19); and (3) analytical calculations performed parallel with 
these model tests using the computer code DIRT (7) with the best possible 
estimate the soil parameters obtained the laboratory tests. Centrifuge tests 
were proposed because them, model scaling relations are adequately satisfied 
covered Refs. and 19. 

The proposal was accepted Chevron Oil Company, and work was initiated. 
The soil involved was clayey silt obtained from the ocean floor off the southern 
California coast. The centrifuge model tests included both monotonic and cyclic 


SIMPLE SHEAR 


Experiment 


sheor stron (%) 


FIG. Barbara Silt: Comparison Prediction and Test Results Simple 
Shear 


loading the model structure, and the results have been reported Refs. 
and 19. This report describes the analytical results and compares them with 
the experimental test results. 


Son 


Soils general undergo both elastic and plastic deformations simultaneously 
upon shearing. order separate the contributions the elastic and plastic 
properties the total deformations, assumed that the elasticity the 
material linear and isotropic, and that nonlinearity and anisotropy result from 
its plasticity. account for the effects unloading-reloading events the 
material behavior, the model combines properties isotropic and kinematic 
plasticity allowing the yield surface translated stress space the 
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stress point and change size simultaneously. Also, order allow for 
the adjustment the plastic hardening rule any kind experimental nonlinear 
hardening law obtained for instance triaxial simple shear soil tests, 
collection nested yield surfaces used. 


stress ratio Tr/ ve 


100 1000 


number of cycles N number of cycies N 


FIG. 3.—Santa Barbara Silt: (a) Derived Cyclic Simple Shear Strain-Controlled Test 
Results; (b) Predicted Cyclic Simple Shear Test Behavior from Analytical Model 
(Figures Curves give Strain percentage) 
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FIG. Three-Dimensional Finite Element Representation Centrifuge 
Model: (a) Side View; (b) Top View Surface 


Since the duration the loading applied the model structure was short, 
was found appropriate assume that deformations the soil foundation 
occurred with volume changes. Therefore, simple version the more 
general formulation reported Ref. was used. For undrained loading 
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FIG. 7.—Vertical Force Versus Vertical Displacement; Comparison Computed 
Results for Three-Dimensional and Two-Dimensional Cases with Experimental Points 
for Centrifuge Model Test 
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FIG. Two-Dimensional Finite Element Representation 
FIG. 6.—Loading System and Notation for Positive Loads and Displacements 
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conditions, the yield surfaces are represented equations the following 
form (14) 


for all which the size the yield surface the deviatoric 
stress subspace; and the center coordinates the yield surface 
the deviatoric stress subspace. Eq. denote the components 
the deviatoric stress tensor. Because the a{” are not necessarily all equal 
zero, the yielding the material anisotropic. Direction therefore 
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FIG. 8.—Computed and Experimental (Test T6) Force Versus Displacement for Three- 
and Two-Dimensional Cases, Inclined and Eccentric Load: (a) Vertical 
Components; (b) Horizontal Components; (c) Moment and Rotation 
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great importance and the physical reference axes (x,y,z) are fixed with respect 
the soil element and specified coincide with the principal axes stress 
the material’s past consolidation history. The initial positions and sizes 
the yield surfaces reflect the past stress-strain history the material element, 
and, particular, their translations are direct expression the material’s 
memory its past loading history associated with the direction 
loading shown Refs. the following the y-axis vertical, 

plastic modulus associated with each the yield surfaces, and 
the associated flow rule used compute the plastic strain rate tensor 
components. The total strain rate set equal the sum the elastic and 
plastic strain rates. All yield surfaces may translated stress space 
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FIG. Force Versus Displacement for Three-Dimensional and Two- 
Dimensional Cases, Inclined Load: (a) Vertical Components; (b) Horizontal Compon- 
ents 


the stress point, and they consecutively touch and push each other but cannot 
intersect. They are allowed change size while they are translated. Their 
associated plastic moduli are also allowed vary, and, general, both 
and are functions the plastic strain history. This further considered 
Refs. and 16. 

Complete specification the model parameters requires the determination 
(1) The initial positions and sizes the yield surfaces together with their 
associated plastic moduli; (2) their sizes plastic moduli changes, both, 
loading proceeds; and (3) the elastic moduli. explained Ref. 14, the 
model parameters required characterize the behavior any given soil can 
derived entirely from the results conventional soil tests. The actual number 
yield surfaces used function the desired degree accuracy 
achieved the material behavior representation. 
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The model parameters required characterize the undrained behavior 
the Santa Barbara silt are described this section. explained Refs. 
and 16, the initial positions and sizes the yield surfaces together with their 
associated plastic moduli, and the elastic shear modulus are determined solely 
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FIG. 10.—Computed Distorted Configuration Soil, Displacements and Velocities; 
Three-Dimensional Test; Inclined and Eccentric Load; 9°; 0.1; 
Side View 


from the experimental behavior the soil measured undrained triaxial soil 
tests. Fig. shows the experimental undrained triaxial stress-strain curves 
obtained for the Santa Barbara silt vertical effective consolidation stress). 
Note that the curves obtained compression and extension are very different. 
Eleven yield surfaces were used fit the experimental stress-strain curves 
closely, and Table lists the parameter values for the Santa Barbara silt, with 
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133.33 The model was thereafter applied predict the results 
simple shear test performed the same soil, and Fig. shows comparison 
stress-strain behavior calculated with the proposed model with that measured 
(9) experimentally undrained simple shear test 1.11 which 
the average shear stress measured the test]. The two sets curves 
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FIG. 11.—Computed Two-Dimensional Solution for Vertical Components 


agree well. These results demonstrate the ability the model account for 
the behavior the Santa Barbara silt upon shearing with rotation the principal 
stress axes. 

the course cyclic loading the yield surfaces all decrease size. 
shown Ref. 16, the change size the yield surfaces cyclic loading 
proceeds determined using the experimental results from undrained cyclic 
strain-controlled simple shear tests. Unfortunately, the only cyclic test results 
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available (9) had been performed under stress-controlled loading conditions. 
The strain-controlled results had therefore derived from the stress-controlled 
results using the procedure explained Ref. 22, and Fig. 3(a) shows the 
corresponding strain-controlled results. The curves shown Fig. 3(a) were 
used determine the cyclic parameters (16) the Santa Barbara silt, and 


FIG. 12.—Computed Two-Dimensional Solutions for Horizontal Components 


Fig. shows the calculated decrease shear stress amplitudes using the 


proposed model. They agree very well with the reconstructed experimental curves 
shown Fig. 3(a). 


The soil model described has been coded and incorporated into finite element 
program (7). With respect undrained loading conditions, should mentioned 
that the finite element program incorporates novel features such penalty 
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function and reduced integration techniques which allow simple and effective 
implementation the incompressibility constraint (10). The use the formulation 
solving boundary value problems illustrated applying analyze the 
interaction the model structure with its soil foundation when subjected 
both monotonic and cyclic forces. The foundation material consists normally 
consolidated deposit Santa Barbara silt, and based the measured average 


FIG. 13.—Computed Two-Dimensional Solutions for Moment and Rotation 


densities the bucket the centrifuge which the model tests were performed, 
equation for the vertical effective stress can written (at prototype 
scale). 


2.50 (0.871 (in tons per square meter) 


which the depth, meters, measured from the surface the soil. 
The (prototype scale) platform has diameter (area 78.54 
and mounted with skirts 0.6 long. Its dead weight such that W/A 
4.3 
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Finite Element Representation.—Fig. shows the three-dimensional finite 
element representation the problem geometry. The complete grid 
the case eccentric loading) spans over 1/2 the bucket geometry, and 
consists binlinear isoparametric brick elements shown Fig. The 
axisymmetric grid (for use the case vertical loading only) spans over 
1/4 the bucket geometry and consists bilinear isoparametric brick elements. 


FIG. Two-Dimensional Distorted Configurations for Inclined and 
Eccentric Loads, Side View, when 10°, ecc/ —0.2 


Three-dimensional calculations are several orders magnitude more expensive 
than their two-dimensional counterparts. Due the necessity working within 
limited computation budget only few three-dimensional calculations were 
performed and most the calculations were performed two dimensions. 
For that purpose the problem geometry was transposed into two dimensions 
assuming that deformation occurred plane strain, and the corresponding 
two-dimensional finite element representation shown Fig. The grid consists 
554 bilinear isoparametric quadrilateral elements. 
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The structure’s foundation represented circular and strip footing 
the three-dimensional and two-dimensional idealizations, respectively. The 
footing consists 0.4-m thick layer elements, and the skirt, where present, 
represented additional layer (row) elements 0.6 deep and 0.4 
wide, attached the footing. The foundation elements are chosen 1,000 
times stiffer than the supporting soil represent the model test conditions. 
Fig. shows the load system and notation. addition the static load 
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FIG. 15.—Computed Two-Dimensional Distorted Configurations for Inclined and 
Eccentric Loads, Side View, when 10° and —0.1 


due the dead weight the structure, the foundation also subjected 
the cyclic inclined eccentric load due the wave (ice) forces. The load 
identified the following means its vertical component acting 
with the eccentricity given rise the moment and its horizontal 
component 

For the case centered vertical loading only the structure 
settles uniformly and denotes its vertical settlement. Under inclined eccentric 
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loading conditions, both, the structure slides and tilts well, shown 
Fig. Its center displaced the amount and the vertical and 
horizontal directions, respectively. The tilt measured the angle rotation 
the slab. 

Analysis Foundation Movements When Subjected Inclined Eccentric Force 
Monotonically Increasing Magnitude.—Figs. and show the computed load 


FIG. 16.—Computed Two-Dimensional Distorted Configurations for Inclined and 
Eccentric Loads, Side View, 10° and ecc/ 0.0 


versus displacement curves for the deformation modes generated various 
inclinations and eccentricities the wave (ice) force under both three-dimen- 
sional and two-dimensional conditions. importance note that for all 
cases, the two-dimensional idealization leads much softer response the 
soil structure interacting system than its three-dimensional (and correct) repre- 
sentation. Also, for comparison, the corresponding centrifuge model test results 
are shown Figs. and and apparent that the three-dimensional analytical 
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predictions agree very well with the measured experimental model test results. 
The experimental results corresponding 10°, 0.0 (Test T2) have 
not been reported Fig. because the poor consolidation state achieved 
this test (19). illustration, Fig. shows the computed three-dimensional 
distorted configuration the soil foundation and the direction its flow. 
apparent comparing Fig. with Fig. the companion paper that 


FIG. Two-Dimensional Distorted Configurations for Inclined and 
Eccentric Loads, Side View, when 10° and 0.1 


the three-dimensional finite element representation does adequately capture the 
failure state displacements the soil foundation. 

Even though appears from Figs. and that only three-dimensional 
calculations are capable representing correctly the three-dimensional model 
test situations, the remaining calculations, due the necessity working within 
limited computation budget, were performed using two-dimensional 
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idealization the problem geometry shown Fig. was felt the 
writers that, though such calculations obviously could not provide exact quantita- 
tive information about the behavior the soil-structure system, they would 
still provide useful answers regarding relative magnitudes loads and displace- 
ments. Figs. 11-13 show the computed load-displacement curves for the deforma- 
tion modes under various inclinations and eccentricities the force and, 
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FIG. Two-Dimensional Distorted Configurations for Inclined and 
Eccentric Loads, Side View, when 1.0° and 0.2 


illustration, Figs. 14-18 show the computed distorted configurations 
the soil foundation and the directions flow failure when 10° and 
for various eccentricities. For inclined and eccentric loading, failure involves 
all three failure modes. failure, very distinct discontinuity occurs the 
velocity magnitude which reveals the existence failure surface. Its trace 
easily obtained from the vectors shown Figs. 14-18 and marked with 
dotted lines these plots. This information can used check the foundation 
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stability using upper bound techniques [see e.g., Ref. (12)], since the latter 
require the determination the failure surface’s trace. the F,, F,, 
space, the various states loading conditions failure describe surface 
which plays the role failure surface for the soil-structure interacting system. 
Its intersection with planes constant shown the interaction diagram 
Fig. for positive eccentricity projections onto the plane 
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FIG. 19.—Failure Surfaces for Soil-Structure System with Positive Eccentricity 
Load 


Analysis Foundation Movements when Subjected Cyclic Inclined and 
Eccentric Force.—Since the behavior the foundation material strongly 
path-dependent, the loading applied onto the structure must follow the real 
wave (ice) loading history closely possible enable predictions 
deformation made. order make the present study specific, the cyclic 
loading applied the structure test (17,19) was selected example 
problem. The computed load-displacement curves for that particular cyclic loading 
case are shown Figs. 20-22 from which the progress the various deformations 
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can traced. importance emphasize that the calculations were 
performed two dimensions and therefore the results not strictly apply 
the real three-dimensional conditions the model test. However, the analytical 
results exhibit the same qualitative behavior observed the test. From 


cyclic 


dead weight 


FIG. 20.—Computed Results Two-Dimensional Cyclic Load Application for Vertical 
Components (Conditions Similiar Test Figs. 14-16 Companion Paper 
but Two Dimensions 


FIG. Results Two-Dimensional Cyclic Load Application for Horizon- 


tal Components (Conditions Similiar Test Figs. 14-16 Companion Paper 
but Two Dimensions 
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Figs. apparent that upon cycling the force the structure accumulates 
permanent vertical and horizontal displacements, and permanent tilt, and that 
the rate accumulation decreases the number cycles loading increases, 
again illustrating the phenomenon (17). apparent comparing 


cyclic +static 


FIG. Results Two-Dimensional Cyclic Load Application for Moment 
and Rotation (Conditions Similiar Test Figs. 14-16 Companion Paper 
but Two Dimensions 


Figs. 20-22 with Figs. the companion paper that the analytical predictions 
exhibit very strong similarity to, and are consistent with the experimental 
results obtained the centrifuge. 


Summary ano 


This investigation concerns offshore gravity structure foundations, their analy- 
sis and design. The results both experimental and analytical studies the 
interaction such structures with their soil foundation when subjected wave 
forces are reported. The experimental tests were performed centrifuge 
100 centrifugal acceleration with 10-cm diam model structure bearing 
normally consolidated deposit Santa Barbara silt. The experimental test 
results therefore apply 10-m diam prototype structure. The model structure 
was subjected both static and cyclic, and inclined and eccentric loading 
conditions, and its resulting movements (sinking, sliding, and tilting) were 
recorded. Although centrifuge tests the behavior gravity platform foundations 
under wave loading conditions have been reported before (23), the former study 
was directed towards the prediction performance specific site with its 
unique soil conditions. the present investigation, point view 
was taken—that validating numerical analysis relatively well-charac- 
terized material. The soil did come from offshore location, but that condition 
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was incidental the tests performed. proved possible carry out the 
investigation the small centrifuge California Institute Technology, 
adequately simulating the response prototype 10-m diam platform soil 
prototype scale properties. Although this not the size the prototype 
production platforms the North Sea, sufficiently large not considered 
model test laboratory scale, and the results are useful their own right. 
The good correspondence experiment and analysis reinforce confidence 
the results both approaches. Since the tests can set and run safely, 
relatively quickly, and reliably the small centrifuge, its use further 
experiments aid assessing the reliability calculational techniques 
encouraging. The analytical predictions for the particular prototype situation 
considered this study were obtained using finite element program which 
incorporates three-dimensional nonlinear soil model capable describing the 
nonlinear, anisotropic, elasto-plastic, path-dependent stress-strain-strength be- 
havior the Santa Barbara silt under both monotonic and cyclic loading conditions. 
The results both two-dimensional (plane strain) and three-dimensional calcula- 
tions are reported and are compared with the appropriate equivalent experimental 
model test results obtained the centrifuge. The three-dimensional analytical 
predictions are found agree very closely with the experimental results and 
the two-dimensional analytical predictions are shown consistent with the 
experimental results. result this study considerable confidence the 
proposed analytical tool has therefore been gained, and felt that should 
now possible use perform reliable analysis actual prototype field 
situations. was shown that the calculations can provide for the structure 
reliable and accurate design information concerning: (1) The overall response 
the foundation (force-displacement and moment-rotation relationships) when 
the structure subjected monotonically increasing inclined and eccentric 
forces various angles inclinations; (2) its various possible failure modes 
(bearing, sliding, tilting) and the trace the failure surface for each case; 
and (3) the overall response the foundation when the structure subjected 
cyclic wave forces. Further, result this study, now clear that 
actual three-dimensional field situations must modeled analyzed such 
accurate experimental analytical predictions are made. 
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INTRODUCTION 


Vibrations produced during construction operations become concern when 
they can felt. They may sufficiently intense annoy the recipients, 
they may have detrimental effects structures equipment. This presentation 
will consider the sources construction vibrations, the transmitting medium, 
and the effects upon the receivers the vibration. Instrumentation and criteria 
will analyzed, well the need for further research. 


Sources Construction 


Construction vibrations are three different types: (1) Transient impact 
vibration; (2) steady-state continuous; and (3) pseudo-steady-state vibrations. 
All are shown Fig. Examples transient construction vibrations are those 
that occur from blasting with explosives, impact pile driving, demolition, and 
wrecking balls. Steady-state vibrations may generated vibratory pile drivers, 
large pumps used jacking underground pipes, and compressors. Pseudo- 
Steady-state vibrations are called because they are random nature 
series impact vibrations that are short enough intervals approach 
essentially steady-state condition. Examples these are jackhammers, pave- 
ment breakers, trucks, bulldozers, cranes, and scrapers. 

the April 14-18, 1980, ASCE Convention and Exposition, held Portland, 
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Ill. 60062. 

Note.—Discussion open until July 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication March 12, 1980. 
This paper part the Journal the Geotechnical Engineering Division, Proceedings 
the American Society Civil Engineers, ©ASCE, Vol. 107, No. GT2, February, 
1981. ISSN 0093-6405 /0002-0167 


167 


FEBRUARY 1981 


Errect Source Enercy 


The magnitude the vibration fixed distance from source function 
the energy the source. This relationship may established increasing 
the energy the source and determining the rate which the vibration increases. 


Types Construction Vibrations 


Transient 


Pseudo 
Steady-State 


Vibration Amplitude Time 


FIG. 1.—Construction Vibration Waveforms 
10.0 


Peak particle velocity 


10,000 100,000 


Energy 
ft-ibs or kg-m 


FIG. 2.—Typical Form: Peak Particle Velocity Versus Energy 


The form typical data (Fig. can approximated straight line 
log-log plot. regression analysis will define the slope this line. Increasing 
the energy can accomplished dropping wrecking ball known weight 
from increasing heights. The energy foot-pounds the product the weight 
and drop height. 
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Although the vibration magnitude may expressed particle displacement, 
velocity, acceleration, has been found most investigators (9,22,31) that 
the peak particle velocity most useful concerns damage and human evaluation 
transient vibration. The relationship between the peak particle velocity and 
energy may expressed 


which peak particle velocity, inches per second; intercept, 
inches per second (value vibration amplitude ft-lb); impact 


10.0 


oO 
> 
= 
x 
a 


1,000 


Distance 


FIG. 3.—Typical Form: Peak Particle Velocity Versus Distance 


energy, foot-pounds; and slope, rate increase. The value has 
been found generally one half. 


TRANSMITTING 


any investigation ground vibration, necessary determine the 
distance propagation equation the rate attenuation the vibration. This 
dissipation energy with distance varies from site site and can determiaed 
recording the vibration simultaneously from several distances, preferably 
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straight line. regression analysis these data can used with propagation 
equation (see Fig. the form 


which peak particle velocity, inches per second; intercept, 
inches per second (value vibration amplitude ft); distance, 
feet; and slope, attenuation rate. 


convenient express the peak particle velocity terms both distance 
and energy single expression. This may done the use scaled-distance 
propagation equation the following form (see Fig. 4): 


which peak particle velocity, inches per second; intercept, 


inches per second [value vibration amplitude D/V 


Pal 
> 
— 
x 
a 


100 


Scaled distance 
| 


FIG. 4.—Typical Form: Peak Particle Velocity Versus Scaled Distance 
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distance, feet; energy foot-pounds for impacts charge weight 


per delay for explosives-pounds; slope attenuation rate; and D/V 
scaled distance. The value this expression generally lies between 
1.0 and 2.0 with relatively common value 1.5. 


Construction 
Typical intensities vibration from the operation construction equipment 


are presented graphically Fig. All the data used prepare Fig. 
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FIG. 5.—Relative Intensities Construction Vibrations 


were obtained field measurements during actual construction operations (29). 
The vibration intensities reported this figure are based data recorded 
the surface the earth residential relatively small commercial 
buildings. The response more massive subsurface structures may 
different, depending frequency input and natural frequency the structure. 

interest note that the intensities involved these construction 
activities are spread over range three orders magnitude. should 
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recognized that this graph represents approximate values, but does illustrate 
the relative intensities the several sources. The lines shown are for particular 
set soil conditions. The locations and slopes lines for other conditions 
may different. 


INSTRUMENTS 


One the early instruments used the 1930s indicate the intensity 
impulse vibrations was the falling-pin seismometer. consisted 1/4-in. (6.4- 
mm) round pins different lengths in.-15 in. (10 cm-38 cm) which were 
placed end level glass plate. Tubes different lengths were bundled 
together that individual pins could topple without knocking down adjacent 
pins. was determined its inventor that transient vibrations which 
the 10-in. pin (25.4 cm) shorter, were capable causing damage 
residential structure. More sophisticated instruments the same era required 
the use tent which modified earthquake seismograph 
was used with light beam drum recorder. 

the 1940s, one the first instruments developed was Leet 
three-component displacement seismograph which utilized system spring- 
suspended masses and light beams. When the light beams were deflected, multiple 
reflections magnified the motion, resulting recording light-sensitive film 
displacement versus time. This system was enclosed aluminum box 
and required 110v 60-cycle current for its operation; weighed Ibs (31.8 
kg). the same period, Knobel accelerograph which recorded particle 
acceleration versus time was developed. This was self-contained battery-operat- 
unit that recorded three mutually perpendicular directions. Ten years later 
(1950), smaller battery-operated displacement recorders were available for field 
use. These generally had magnification and weighed approximately 
Ibs (15.9 kg). 

Vibration recording instruments common use since 1960 utilize coils 
suspended weights soft springs with magnetic core attached the frame 
the instrument. The voltage generated when the frame vibrates function 
the particle velocity the vibration. means light beam galvanometers, 
recording particle velocity versus time obtained; several sensitivity levels 
can selected. The more sophisticated instruments automatically trigger and 
identify the time and date the event. relatively inexpensive strip chart 
recorder which gives bar graph the maximum vibration also available 


for long-term monitoring, are magnetic tape analog equipment and digital 
recorders. 


Crrreria For Limitine Viprations 


Vibrations from construction operations are important primarily because 
their potential cause complaints damage variety structures, 
both, above and below ground. They may also cause damage building contents 
such china collections, exhibits museums, more commonly, sets 
glasses. Vibrations have resulted complaints that they are disturbing 


people the enjoyment their homes interference with the performance 
their jobs. 
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The most frequent complaints are from the owners structures (particularly 
residential) who claim that damage has been done their interior walls 
exterior walls and foundations, both. probably for this reason that most 
attention has been focused residential-type construction. Investigators 
this continent and abroad have criteria for limiting such vibration. 
Some are simple, some are quite complex. The United States Bureau Mines 
(USBM) (22) and Edwards and Northwood (9) recommended safe levels 
2.0 ips (50.8 mm/s) for residential structures (see Fig. 6). One investigator 
(4) classifies buildings four categories based their age and condition, and 
recommends criteria (Table that range from 0.5 ips (12 
mm/s), the former for old residential structures very poor condition, the 


SERIOUS 
CRACKING 


MAJOR DAMAGE 
(FALL OF PLASTER, 
SERIOUS CRACKING) 


CRACKING 
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OF OLD CRACKS) 
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CAUTION 


Particle 
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DAMAGE 
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LANGEFORS EDWARDS BUMINES 
(SWEDEN) (CANADA) (U.S.A.) 
(17) (9) (22) 


FIG. 6.—Comparison Damage Criteria: Residential-Type Structures 


latter for structures substantial construction. Another investigator (6) recom- 
mends criteria that take into account whether the building ancient modern, 
specific construction details, and degree distress for both ancient and modern 
buildings, type building, type soil, the cumulative effect vibration, 
and factor.’’ The latter ranges from light industrial sensitive, 
ancient, historical buildings, scale 1-10. 

The new Swiss standard for vibration (32) differentiates between machines, 
traffic, and blasting, categorizing all buildings four classes. For each class 
building, permissible levels vary with frequency, shown Table The 
Swiss criteria are more restrictive than the for Langefors (17), Edwards 
and Northwood (9), and the United States Bureau Mines (22). 

relate construction vibrations damage criteria, Fig. presented. Fig. 
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TABLE 1.—Criteria Recommended Chae (4) 


Peak particle 
velocity, 
millimeters 
per second 


(3) 


substantial construction. 

new residential structures sound condition. 

“Relatively old residential structures poor condition. 

“Old residential structures very poor condition. 

Note: structure subjected repeated blasting, blasting done without 
instrumentation, lower category one. 


TABLE 2.—New Swiss Standard for Vibrations Buildings (32) 


Peak particle Peak particle 
Range velocity, velocity, 
Building Vibration frequency, millimeters inches 
class source hertz per second per second 
(1) (2) (3) (4) (5) 
Machines, traffic (0.5) 
(0.5-0.7) 
Blasting (1.2) 


(1.2-1.6) 
Machines, traffic (0.3) 


(0.3-0.5) 

Blasting (0.7) 

Machines, traffic (0.2) 
(0.2-0.3) 

Blasting (0.5) 
(0.5-0.7) 

Machines, traffic (0.12) 

(0.12-0.2) 

Blasting (0.3) 


(0.3-0.5) 


steel reinforced concrete, like factories, retaining walls, bridges, steel 
towers, open channels; underground chambers and tunnels with and without concrete 
alignment. 

with foundation walls and floors concrete, walls concrete masonry; 
stone masonry retaining walls; underground chambers and tunnels with masonry alignments; 
conduits loose material. 

mentioned previously but with wooden ceilings and walls masonry. 

“Construction very sensitive vibrations; objects historic interest. 


Class 
(1) 
100 
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the same Fig. previously presented, but has superimposed upon 
the levels vibration intensities which generally are considered investigators 
the United States thresholds possible damage for both commercial 
buildings and for residential-type structures. These are 4.0 ips for commercial 
buildings and 2.0 ips for residential structures. should recognized that 
some the sources are transient (impulsive) vibrations, others represent 
steady-state, and others are pseudo-steady-state vibrations. For purposes 
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FIG. 7.—Construction Vibration Intensities Related Damage Criteria 


this presentation, they have been lumped together. emphasized that these 
are which generally must exceeded before damage 
likely occur. 

generally acknowledged that the ratio the frequency the vibrations 
transmitted the natural frequency the structure important parameter. 
The closer this ratio unity, the greater the response the structure, 
therefore the lower should the criteria. Since residences have natural 
frequencies the range Hz-10 Hz, lower criteria have been suggested 
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for lower frequency ground vibration than for higher frequency vibration. 

Several investigators (7,20,21) have been and are working the response 
spectrum. this method, the amplitude the vibrations induced structure 
function the natural frequency and damping characteristics the structure 
and the frequency composition and magnitude the ground motion. 

Published criteria for underground utilities are not available. Such utilities 
include tile concrete drainpipe, sewer and water mains, and pipelines. 
our experience, particle velocities ips have not damaged pipes and mains. 
Further, high pressure pipelines have withstood ips without experiencing 
any distress, apparent from dynamic strain gage measurements. 
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(28 day) 
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@ 
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Curing time days 
FIG. 8.—Strength Concrete Versus Curing Time 


Another problem that frequently arises construction jobs the safe level 
transient vibration the vicinity young concrete. The primary sources 
such are pile driving and excavation blasts. With young concrete, are 
dealing with material which increases strength with time, thus, must 
the criteria. The writer (27) considers the concrete essentially full 
strength after days, Ontario Hydro (unpublished) days. decision must 
then made the particle velocity that full-strength concrete can withstand. 
The writer (27) says ips, Atkins and Dixon (1) ips, Ontario Hydro (unpublished) 
ips. Because stress level proportional particle velocity, the criteria can 
established when known how the concrete attains strength with time 
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(see Fig. 8), e.g., when the concrete attains 10% its strength, the permissible 
level for vibration would 0.5 ips, 0.4 ips, 0.3 ips. When attains 50% 
strength, the permissible level would 2.5 ips, 2.0 ips, 1.5 ips. 

The effects transient vibrations upon soils another area that concern 
engineers. Some work has been done the effect steady-state vibrations 
soils (24). Little, any, quantitative information available the effects 
short-duration transient vibrations soils. best, the literature (5) states 
that cohesive soils are least affected, e.g., compacted, vibration, while 
cohesionless soils are most sensitive vibration. 
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FIG. 9.—Comparison Human Response Steady State (23) and Transient Vibra- 
tions (28) 


Information available the effect cyclic loading various materials. 
These have been experimentally determined for complete reversals stress 
several magnitudes and for partial reversals stress. The studies have involved 
repetitious steady-state reversals. There information repeated transient 
loading. The difference that transient, the peak stress occurs but once 
twice with the remaining cycles being reduced magnitudes. far 
known, there virtually published information fatigue studies any 
kind plaster masonry materials used building construction. 

investigation into the fatigue plain concrete (15) reveals that the fatigue 
strength for life 10,000,000 cycles roughly about 50% static strength 
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for large-amplitude stress cycles, and approximately 85% for smaller cyclic stress 
amplitudes. Reinforced concrete exhibits practical fatigue limit the vicinity 
1,000,000 cycles and endurance limit excess ksi (1,409 kg/sq 


cm). 


Human Response 


Regardless the criteria that apply transient vibration, rare 
that the criteria can reached before complaints claims are filed. The 
human body excellent detector vibration, but poor measuring device. 
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FIG. 10.—Construction Vibration Intensities Related Human Evaluation Vibra- 
tions (23) 


Considerable work has been done the human evaluation steady-state 
vibrations (14,23), but very little has been done transient vibrations. The 
only study found transient vibrations (28) was concerned with impact vibrations 
which simulated footfalls while walking across floor. The waveform 
sufficiently similar those resulting from impact pile driving from single-hole 
blasts that the criteria would apply. Multidelay blasts not have waveform 
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that similar the footfall vibrations. The human evaluation the transient 
vibration study compared those ratings reported from steady-state vibration 
studies shown Fig. 

relate construction vibrations human response steady-state vibrations, 
Fig. presented. Fig. the same Fig. previously presented but 
has superimposed upon the normal subjective evaluation the vibration 
intensity. apparent comparison Figs. and that the vibration 
intensities which are subjectively classified 
are well below the intensities required cause damage. Levels the 
range may result complaints claims damage. 


considered that the instrumentation needed measure and document 
vibration during construction operations presently available. Although expen- 
sive, competition will bring the costs down. 

The development yet not overly conservative criteria needed. 
This specifically applies residences, buildings (both low-rise and high-rise), 
and other structures, well underground utilities. Those currently use 
are either too simple too complex. The criteria should based available 
and future controlled tests and not arbitrary evaluations. 

Studies evaluate the criteria used for young concrete can made 
the field the laboratory simulating transient vibrations hydraulic 
shaker. The effects transient vibration soils can studied the field 
with laboratory tests providing support information. 

The effect fatigue common building materials, particularly plaster and 
masonry, should investigated. This should done for both continuous stress 
reversals, well transients. acknowledged that, when dealing with 
transients, the problem more complex. When dealing with continuous reversals, 
only the frequency, amplitude, and time are parameters. With transients, there 
are the additional parameters damping, duration individual pulses, rate 
occurrence, and complex waveforms. 

Human response transient vibration being studied research program 
entitled Response Blast and Noise and Ground and Structure 
presently underway contractor for the United States Bureau 
Mines. This 2-yr program which consists attitude surveys and field 
and laboratory tests. should provide meaningful results the construction 
industry, even though the blasts investigated the USBM program will 
surface coal mines and quarries. 


Summary 


The generation and propagation vibrations due construction have been 
considered. The relative intensity the vibration produced many construction 
operations and equipment are identified, based more than quarter century 
field tests. Instruments suitable for the measurement these vibrations 
have been described. The large differences current damage criteria are revealed. 
Research needed the effects vibration, both transient and steady-state, 
young concrete, soils, and underground utilities. Fatigue data building 
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materials for both transient and steady-state vibrations are needed. Additional 
research needed the human response transient vibrations. 

The state the art has come long way the last yr. With sufficient 
funding, all the required research can done; priorities will established 
the need. Hopefully, will not take the next for the required research. 


Atkins, P., and Dixon, E., Structures and Construction 
Vibrations Concrete Structures, SP-60, American Concrete Institute, 1979, pp. 
213-247. 

Attewell, B., and Farmer, W., Ground Vibrations from Pile 
Ground Engineering, Vol. No. July, 1973. 

Brown, M., Vibrations Caused Construction Equipment 
and Report No. RR172, Department Highways, Ontario, Canada, Apr., 
1971. 

Chae, S., Excavation Blasts Prevent Civil Engineering, 
ASCE, Vol. 48, No. Apr., 1978, pp. 77-79. 

Chellis, B., McGraw-Hill Book, Co., Inc., New York, N.Y., 
1961. 

Crockett, A., Vibrations and Structural Fatigue,’’ Recent Developments 
the Design and Construction Piles, Institute Civil Engineering, London, England, 
1979. 

Dowding, H., Response Blasting Proceeding 
the 14th Annual Meeting, Society Engineering Sciences, Lehigh University, Nov., 
1977, pp. 1085-1098. 

Dym, L., Ground Vibration,’’ Sound and Vibration, Apr., 1976, 
pp. 32-34. 

Edwards, T., and Northwood, D., Studies the Effects 
Blasting The Engineer, Vol. 210, Sept., 1960, pp. 538-546. 

Ferahian, H., and Hurst, D., and Possible Building Damage due 
Operation Construction Proceedings the 1968 Public Works 
Congress and Equipment Show, Oct., 1968, pp. 144-155. 

Ferahian, H., Caused Pile Bibliography No. 35, National 
Research Council Canada, Division Building Research, Ottawa, Canada. 

Ferahian, H., and Ward, S., Environment Laboratory 
Technical Paper No. 329, Division Building Research, National Research Council, 
Ottawa, Canada, Oct., 1970. 

Forssblad, L., Vibrations and the Damage They Cause, the Effects 
Vibratory Compactors,’’ Swedish Geotechnical Institute, May, 1973. 

Goldman, A., Review Subjective Responses Vibratory Motion 
the Human Body Frequency Range Report NM-004-00, Naval Medical 
Research Institute, Washington, D.C., 1948. 

Hanson, M., Ballinger, A., American Concrete Institute Committee 215, 
erations for Design Concrete Structures Subjected Fatigue Loading,’’ Journal 
the American Concrete Institute, Vol. 71, No. Mar., 1974, pp. 97-121. 

Heckman, S., and Hagerty, J., Associated with Pile 
Journal the Construction Division, ASCE, Vol. 104, No. CO4, Proc. Paper 14205, 
Dec., 1978, pp. 385-394. 

Langefors, U., B., and Westerberg, H., Vibrations 
Water Power, Feb., 1958. 

Liu, K., and Kinner, B., Yegian, K., Sound and 
Vibration, Oct., 1974, pp. 26-32. 

Luna, A., Vibration Due Pile Foundation Facts, Vol. 
No. Raymond International, Houston, Tex., 1967. 

Medearis, K., Development Rational Damage Criteria for Low-Rise Structures 
Subjected Blasting Proceedings the 18th U.S. Symposium Rock 
Mechanics, June, 1977, pp. 1A2-1-1A2-6. 


3 


GT2 CONSTRUCTION VIBRATIONS 181 


21. Naik, R., Damage Low-Rise Buildings Due Ground Vibrations 
Created Vibrations Concrete Structures, SP-60, American Concrete 
Institute, 1979, pp. 249-264. 

22. Nicholls, R., Johnson, F., and Duvall, Vibrations and their 
Effects BuMines Bulletin 656, 1971. 

23. Reiher, H., and Meister, J., Effects Vibration (In German) 
Forschung auf dem Gebeite der Ingenieurwesens, Vol. No. 1931, 381. 

24. Richart, E., Hall, R., and Woods, D., Soils and 
Prentice-Hall, Inc., Englewood Cliffs, N.J., 1969. 

25. Steffens, Aspects Structural Proceedings the Symposium 
Vibration Civil Engineering, Imperial College Science and Technology, Apr., 
1965. 

26. Sutherland, B., Study Vibrations Produced Structures Heavy 
Proceedings the Thirteenth Annual Meeting the Highway Research Board, 
1950, pp. 406-419. 

27. Wiss, F., Effects Pile Driving Highway Research Board 
Record 155, 1967, pp. 14-20. 

28. Wiss, F., and Parmelee, A., Perception Transient 
Journal the Structural Division, ASCE, Vol. 100, No. ST4, Proc. Paper 10495, 
Apr., 1974, pp. 773-787. 

29. Wiss, F., During Construction Journal the Construction 
Division, ASCE, Vol. 100, No. CO3, Proc. Paper 10798, Sept., 1974, pp. 239-246. 

30. Wiss, F., Much Vibration,’’ World Construction, Oct., 1974, pp. 24-26. 

31. Wiss, F., Study Damage Residential Structure from Blast 
Special Publication, ASCE, 1974. 

32. Vibration Swiss Association Standardization, See- 
feldstrasse 8008, Zurich, Switzerland. 


The following symbols are used this paper: 
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coefficients; and 


exponents. 
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RESPONSE BURIED STRUCTURES 
TRAVELING WAVES 


The response underground structures, mainly pipelines, receiving more 
and more attention because the impact such lifelines water mains and 
sewer lines upon the health and safety people served and also because 
the hazardous consequences the failure gas oil lines. Procedures are 
available for analyzing these structures subjected earthquake loadings (1,2,3, 
6,7,11,13). However, all these procedures were developed the assumption 
that soil-structure interaction effects the response structures are negligible, 
i.e., structures move with the surrounding soils during earthquakes. 

reasonable assume that the soil-structure interaction effects are negligible 
for small and flexible structures. For large structures, such subways, the 
potential errors introduced ignoring the interaction effects can certainly not 
overlooked. the object this study evaluate the interaction effects 
the response large structures traveling seismic waves. new approach 
introduced for computing the response these structures using the finite 
element method. analyses using this new approach the soil-structure systems 
are represented three-dimensional finite elements; however, using special 
type wave-transmitting element, three-dimensional system can reduced 
plane finite elements with three degrees-of-freedom each node while 
the traveling effects waves are faithfully retained. 
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Fig. shows typical interaction system interest. The system 
represented three-dimensional finite element model. For particular cases 
which the geometry and material properties remain unchanged along the 
longitudinal axis the structure, the system can modeled into identical 
planes finite elements equal spacings. All these planes are parallel 
each other and perpendicular the longitudinal axis the structure shown 
the figure. 

Consider brick element bound two parallel planes and shown 


FIG. 1.—Finite Element Representation Soil-Structure Interaction System Typical 
Buried Long Structure: (a) Three-Dimensional Model; (b) Wave-Transmitting Element 


and the corresponding nodes plane and For simplicity only 
three translational degrees-of-freedom each will considered. 
stiffness matrix [K] can formed for this element and the forces 
acting these nodes can related the nodal displacements 


which and vectors containing forces acting nodes planes 
and respectively; similarly, and nodal displacements these nodes; 
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this element alone, are 


For waves traveling the longitudinal direction constant velocity, the 
following relationship exists: 


which the wave number; and the thickness the element 
the spacing between two consecutive planes. 

The wave number, related the velocity, and frequency, the 
wave the following expression: 


should noted that the wave number used here global wave number 
refering the frequency and the velocity the wave propagating through 
the entire profile and does not relate the wave velocity the material making 
the element. However, readers are reminded that assumption that the 
geometry and the material properties remain unchanged along the longitudinal 
direction the structure has been made for Eq. valid. This aspect 
further considered the end this paper. 
Substitution Eq. Eq. yields: 


Similarly, the additional forces acting nodes plane due the deformations 
the neighboring element bound planes and [see Fig. are 


The total forces acting nodes plane can thus obtained combining 
these two equations, giving 


The preceding expression relates the forces and the displacements the 
four nodes single plane; therefore, definition, the stiffness matrix for 
the wave-transmitting element formed these four nodes 


The stiffness matrices for all the elements the same plane can thus formed 
and the global stiffness matrix can easily obtained assembling the element 
stiffness matrices. The three-dimensional system therefore reduced plane 
finite elements. 

can noted that although the submatrix [K,,] symmetric matrix, 
the submatrices and are not; therefore the element stiffness matrix 
[L] not longer symmetric matrix. Furthermore, because the presence 
the imaginary terms, [L] not even real matrix, the one commonly 
used plane strain analyses. matter fact, because the two elements 
the two sides the plane are identical 
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which the terms parentheses are the corresponding terms the submatrices 
and [K,.], respectively. The element stiffness matrices, and 
therefore the global stiffness matrix, are thus Hermitian matrices material 
damping ignored. 

The foregoing procedure straightforward and easy understand. However, 
difficult form the submatrices [K,,] and for arbitrary brick 
element without first forming the full element stiffness matrix [K] which 
matrix. The formation full stiffness matrices for all the elements 
very costly. simpler method for forming the stiffness matrices for wave- 
transmitting elements given the next few paragraphs. 

standard practice forming the stiffness matrix for quadrilateral element 
transform the x-y system new coordinate system that the coordinates 
arbitrary point the element can expressed terms the coordinates 
the four nodes the corners, thus 


4 a 

i=l i=l 


coordinates the ith node the element. Similarly, the displacements can 
also expressed terms nodal displacements 


For waves traveling the longitudinal direction, the direction, the displacement 
field for the entire system can expressed 


i=l 


Therefore 
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From calculus, the following relationship exists: 
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Differentiation Eq. yields 


(u, orw) [T] 


T T 
T;, Tx 


Strains are related displacements 


— +t) 
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Using the relationships Eqs. and 21, one can obtain the following expression: 


{e} 


inwhich 


For linear elastic isotropic material the stress-strain relationship given 


Hooke’s 


Law 


which 


€y 
(24) 
{u} 
T\; T 3, 7 | 
(26) 
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where the Young’s modulus; and the Poisson’s ratio. 
Using the virtual work approach the expression for the stiffness matrix 


The integration usually carried out using Gauss quadrature formulas which 
turn the integration into the summation finite number terms. should 
carried out through the entire volume the element. However, possible 
form the matrix such way, see Eq. that 


[B(s, 


the stiffness matrix for element like the one shown Fig. can then 
obtained follows: 


2id area 

which [B*] the transpose conjugate (or called tranjugate) [B*]. 


these slices are sufficiently thin, implying continuous system, this equation 
reduces 


which [L*] dA. Because only the four nodes 
plane are involved the formulation, the matrix [L*] can imagined 
the stiffness matrix for plane element bound these four nodes while 
waves are free pass through. The formulation element stiffness matrices 
using this method much easier than the formulation using Eq. because 
the smaller number degrees-of-freedom involved. also more rigorous 
because Eq. exact representation the variations displacements 
the longitudinal direction once the assumption that the geometry and the 
material properties remain unchanged made. forming the stiffness matrix 
for arbitrary three-dimensional element, usually assumed that displace- 
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ments vary linearly with space. Therefore, strictly speaking, Eq. not exact. 
However, the thickness element, approaches, the two methods approach 
the same results. 

should noted that, because the presence the imaginary terms 
the third row the matrix [P], the matrix [B*] complex matrix and 
matrix [L] longer real and symmetric. has been shown previously, 
the stiffness matrix [L] now Hermitian matrix material damping ignored. 
This becomes clear after the matrix operation Appendix followed. 

Material damping can easily incorporated analyses using the complex 
response method (5). The stiffness matrix for element for which the damping 


[L] 


such cases the matrix [L] longer Hermitian. 
Motion 


After the stiffness matrices for all the elements plane have been formed, 
global stiffness matrix can obtained assembling all these element stiffness 
matrices. The equation motion for slice elements thus 


which [M] the mass matrix; [K] the stiffness matrix; and {p} 
external forces. The mass matrix [M] the same the one usually seen 
plane strain analyses and can formed using the lumped-mass method, 
the consistent-mass method, the combination the two. The external forces 
are zero free nodes and nonzero nodes the base the finite element 
model. Given the displacements nodes the base, the motions the free 
nodes can obtained solving the equation 


which reduces 


The subscript the foregoing expression refers free nodes and the 
nodes the base. This equation commonly seen and techniques are available 
for solving for response (5). Because the presence complex numbers, 
transient response will best solved the frequency domain using the Fourier 
techniques decomposing the input motions into serious harmonic motions 


j=l 
j=l 


The same treatment can applied the response 


j=l 
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j=l 


the Fourier amplitude input displacement; circular frequency; time; 
and the Fourier amplitude response. Therefore, for specific frequency 


This set simultaneous linear equations can solved using either iterative 
elimination schemes for {U,}. The total response can then obtained using 
and 39. 


VERIFICATION THEORY 
computer code STABLE (Simplified Three-Dimensional Analyses Buried 


Long structures during Earthquakes) has been developed carry out the 
computations outlined the previous sections. verify the formulations, the 
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FIG. Problem: Propagation Rayleigh Waves 
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response three layer system, Fig. was computed and compared with 
that obtained the previous studies. The system was first analyzed Stoneley 
(9), later Stoneley and Hochstrasser (10), and again the second writer 
(4). 

The results obtained using the code STABLE are compared with the results 
published Stoneley Table this particular case analyzed, harmonic 
wave travels phase velocity 12,441 fps (3.79 km/s) and frequency 
2.48 Hz. can noted, the agreement between the two sets results 
excellent. Since Stoneley’s results are closed-form solutions, the comparison 
provides vital test the proposed theory. 


TABLE 1.—Computed Response for Problem 


From Stoneley (9) This Study 


Depth, feet Horizontal Vertical Horizontal Vertical 
(1) (2) (3) (4) (5) 


Solutions for more complicated problems, unfortunately, are not available. 
Therefore, although desirable, not possible test the code for problems 
such soil-structure interaction problems for which the code was developed. 


Demonstration 


demonstrate the potential applications the proposed procedure the 
response 28-ft 20-ft (8.5-m tunnel was analyzed. The lining 
the tunnel made 2-ft (0.6-m) thick reinforced concrete for which the 
Young’s modulus 5.2 10° psi (0.366 10° kg/cm’) and the Poisson’s 
ratio 0.24. The finite element model used the analyses shown Fig. 
Also shown the figure are the soil profile and soil properties. 

The time history the input motion the base the finite element model 
shown Fig. artificial acceleration time history generated 
accordance with Regulatory Guide 1.60 spectra (8) give broad frequency 
content the frequency range interest. Although possible specify 
input motions all the three directions, for this particular case, the motion 
assumed horizontal and travels the longitudinal direction the 
tunnel. The phase velocity the waves assumed 2,000 fps (610 m/s). 

fully investigate the effects wave travelling the response the 
soil-structure system, the following cases have been analyzed: 


—0.8049 1.0000 —0.8043 0.9988 
200 —0.5800 1.0313 —0.5793 1.0304 
400 —0.4034 1.0254 —0.4028 1.0245 
600 0.9895 —0.2728 0.9888 
800 —0.1890 0.9312 —0.1887 0.9307 
900 0.9097 —0.1436 0.9094 
1,000 —0.1089 0.8848 —0.1088 0.8846 
1,100 —0.0841 0.8572 —0.0841 0.8572 


GT2 BURIED STRUCTURES 193 


Free field soils only, assuming the foundation rock perfectly rigid. 

Soil-tunnel interaction system, assuming the foundation rock perfectly 
rigid. 

Free field soils only, assuming the foundation rock flexible and waves 
are traveling velocity 2,000 fps (610 m/s). 

Soil-tunnel interaction system, assuming the foundation rock flexible 
and waves are traveling velocity 2,000 fps (610 m/s). 


The assumption that the foundation rock perfectly rigid implies infinite 
speed wave traveling. such cases, the analyses are essentially the same 
the conventional two-dimensional plane strain analyses commonly performed. 
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Therefore, the computer code FLUSH (5) was used compute the response 
for these cases. 

The distributions peak accelerations the free field and the transfer 
functions, defined the ratios between the computed response and the input 
motion function frequency, the motions the surface are shown 
Fig. The effects wave traveling can readily seen. the two cases 
shown, vertical response was induced. 


Horizontal Acceleration at Surface 
Horizontal Acceleration at Base 


Frequency - Hz 


Peak Acceleration - 9 
0.5 1.0 


FIG. 5.—Free Field Response 


Fig. shows the transfer functions the horizontal components the response 
obtained depth (6.1 m), which the top the tunnel, for 
the first four cases analyzed. The upper figure for the two cases which 
the foundation rock assumed flexible (traveling wave assumption) and 
the lower figure for the two cases which the foundation rock assumed 
rigid (vertically propagating wave assumption). can noted, the 
soil-tunnel interaction effects the response are not significant. However, 
the wave-traveling effects, can evaluated comparing the upper and 


C = inf. | 
An 
/ 
/ 
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the lower figures, the response can not ignored. general, can 
concluded that, traveling waves are assumed, peaks transfer functions 
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CENTERLINE OF TUNNEL 


Input MoTIoN AT BASE 


FREE-FIELD 


Frequency - Hz 
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FIG. 7.—Transfer Functions Shear Strains Depth (6.1 


shift toward higher frequencies and the amplitudes the peaks become less. 
This conclusion agrees with previous observations (12). For the cases studied, 
the vertical response caused soil-tunnel interaction quite small. 
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The transfer functions the shear strains, y,, and computed the 
upper corner the tunnel are compared with those induced the same depth 
the free field Fig. Contradictory the conclusion just made that soil-tunnel 
interaction effects are insignificant the computed response, the presence 
tunnel drastically reduces shear strains; therefore, will overly conservative 
design the tunnel lining based free field soil strains. 

The fact that different conclusions were reached studying the motions 
and the strains seems surprising but totally reasonable. For example, replacing 
thin soft layer extra stiff layer horizontally layered system subjected 
horizontal shaking will hardly make any differences the computed motions. 
But the shear strains induced this extra stiff layer will much less than 


Peak STRESS KSF 
5.0 7.5 


=x 
a 
WW 


FIG. 8.—Shear Strains Free Field Versus Depth 


those induced the original soft layer. This because the fact that strains 
are differentiation displacements with respect distance while accelerations 
are double differentiation displacements with respect time. Strains and 
acceleration are usually not directly related except for particular cases. the 
example mentioned, although the computed motions are almost the same for 
the two systems, they are not identical. Because the stress equilibrium 
the interfaces between adjacent layers, the strains induced the extra stiff 
layer and the original soft layer will inversely proportional the shear 
moduli the two layers. Therefore will wrong reach conclusions 
the significance soil-structure interaction effects comparing motions 
obtained strains, instead motions, are the design parameters interest. 
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Unfortunately, this important fact has not been well recognized the past. 
Another observation which can made studying Fig. that the traveling 
wave assumption yields significant amount which not exist vertically 
propagating waves are assumed. Accordingly, the assumption vertically 
propagating waves may not adequate this case. 
Fig. shows the distributions the computed shear strains the free field 
versus depth. Pronounced wave-traveling effects can easily noted. The peak 


TABLE 2.—Comparison Strains Induced Free Field and Structural Response 


SHEAR STRAINS, percentage 


Free Field Corner Tunnel 


Depth, feet Component 
(1) (2) 
0.067 
0.106 


Peak ACCELERATIONS - G 


C= 2000 rps = INFINITY 


Horiz, 
VerT. 
AXIAL 


FIG. 9.—Peak Structural Response 


shear strains obtained the corners the tunnel are summarized Table 
Also shown the table are the corresponding free field strains. general 
the strains the corners the tunnel are less than one fifth those 
the free field except for which the free field values are zero. 

The peak stresses computed the corners the tunnel are listed Fig. 
interesting see that the traveling wave assumption also yields high 
normal axial stresses result the bending the tunnel waves 
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(4) (5) (6) 
0.044 0.013 0.007 
0.002 
0.067 0.004 
0.064 0.014 0.007 
0.001 
0.060 0.003 
PEAK STRESSES PSF 
0.08 33.800 N/A 
Vert. 0.03 0.04 16 } 
12.370 
14' 14' Txy 21,900 40.970 
x 
z 
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travel along its longitudinal axis. For vertically propagating waves, some amount 
axial stresses are caused Poisson effects because the plane strain 
assumption. Their magnitudes are not readily available but are limited 
which gives 7,340 psf (351,586 Pa) the upper corners and 6,970 psf 
(333,863 Pa) the lower corners. These amounts are significantly lower than 
those obtained for traveling waves. 

expected, the presence tunnel produces some amount vertical motions 
and the traveling waves also causes some longitudinal movements. However, 
these induced motions are relatively small compared the horizontal motions 
because this case for SH-wave propagation. 


The proposed approach subjected couple limitations. First all, 
the geometry and the material properties have remain the same along the 
longitudinal axis structure order for Eq. valid. This rarely 
the case. Therefore, most cases, representative soil profiles will have 
used. This line with almost all the analyses involving soils, even static 
problems, and should not present any difficulties for practical purposes. Secondly, 
the analyses using this approach are essentially linear analyses because responses 
are solved the frequency domain instead the time domain. The nonlinearity 
material properties will have incorporated the analyses using the 
equivalent linear method conjunction with iterative scheme. 

The proposed approach can expanded analyze the response soil-struc- 
ture interaction systems surface waves using frequency-dependent wave 
numbers which can determined solving the eigenvalues the free field 
system obtain the natural frequencies the system. Methods are available 
for such task (4,14). The approach can also used, with minor modifications, 
for other types loadings such moving trains rails, landing aircrafts 
etc. The structures not have underground. Surface structures 
such long dikes, retaining walls can also analyzed equally well. 


The approach presented herein provides useful and practical tool for studying 
the respone soil-structure interaction systems traveling waves. Based 
the results obtained, the following concluded: 


For large structures, although the interaction effects the computed 
motions are small, the presence the structures significantly reduces the 
computed strains because the rigidities the structures provide much resistance 
deformations. Therefore, may overly conservative design such 
structures based the strains developed the free field. 

The assumption vertically propagating waves are not adequate for 
predicting stresses and strains induced underground structures. Proper consid- 
erations should given the traveling effects waves the design. 


The multiplication matrices Eq. carried out herein. Let the matrix 
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[B*] separated into the real and the parts that 

then the multiplication can carried out follows 


Because [H] real and symmetric matrix, the real part the result shown 


the right hand side this expression will symmetric matrix. Furthermore, 
because the fact that 


the imaginary part the result will antisymmetric matrix. Therefore, 
the product [B*] Hermitian matrix stated previously. 
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The following symbols are used this paper: 


area; 

phase velocity waves; 
Young’s modulus; 
forces, frequency; 
thickness slices; 
Jacobian; 

stiffness matrices; 
stiffness matrices; 

mass matrices; 

forces; 

transformed coordinates; 
absolute displacements the and directions; 
volume; 

original coordinates; 
damping ratio; 

strains; 

wave number; 

Poisson’s ratio; and 
stresses. 
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LATERAL PRESSURE PEAT SOILS 


INTRODUCTION 


The coefficient lateral earth pressure under conditions lateral 
deformation, has been interest geotechnical engineers for many years 
for its role reconstructing the geotechnical history soil deposits and for 
its being essential parameter the design analysis many conventional 
problems. The number situations requiring knowledge has increased, 
especially during the last decade two, result the increase types 
structures and situations which now require great deal more rigour and 
thoroughness their analysis. Although not readily amenable mathematical 
formulation terms other soil parameters, complex function 
soil microstructure (fabric), strength, and stress history. These parameters are, 
general, interdependent certain extent. The range values K,, and 
its relationship number basic soil parameters, have been reported for 
sand and clay soils (see Abdelhamid and Krizek (1) for comprehensive list 
references); however, there are very few reports for peat soils (2). 
Peat, encountered extensively many parts the world, fragmented organic 
material derived from vegetation which has been chemically changed and 
fossilized. formed wetlands under appropriate climatic and topographic 
conditions. Today, some the marginal sites are old wetlands that have been 
drained. Because development some parts the world where peat deposits 
are extensive, methods improving engineering properties peats situ 
and thus utilizing them foundation materials, are becoming more and more 
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common (12). prime importance for the success such procedures 
understand and quantify various engineering properties peat soils. 
attempt made this direction, this study presents the results experimental 
program directed towards evaluating and its dependence stress history 
and other parameters peat soils. 


Peats 


Fiber content appears major compositional factor determining the 
way which peat soils behave. For this program, three types peat with 
fiber contents ranging from about 20%-60% were investigated. The sampling 
sites the peats used this study are located the following places 
Wisconsin: Fon Lac county, the city Portage, and the city Middleton 
(5), and the peat samples are designated after these locations. Fond Lac 
Peat amorphous granular peat which the soil particles are mainly 
colloidal size and the other two peats are fibrous peats which have essentially 
open structure with interstices filled with secondary structural arrangement 


TABLE 1.—Average Properties Peat Samples 


Unit 
weight, 
in kilo- 
newtons 


Classification 


Vane 
shear 


MacFarlane 


per strength, and 
Source of cubic in kilo- Sensi- Radforth 
samples meter pascals tivity (Ref. 13) 


(1) 
Fond du Lac 


(3) (8) (9) (11) 


Amorphous 


County, Wisc. granular 
Portage, Wisc. Fibrous 
Middleton, Wisc. Fibrous 


nonwoody fine fibrous material (13). all the sites, the peat deposit 
occurred the ground surface, and sampling was achieved using 76-mm 
Shelby tube (thin-wall sampling tube with very sharp leading edge) (5). 
portion the sample each tube was used determine the characteristics 
the peat. The average properties obtained for the three types peat used 
this study are given Table These properties were determined largely 
accordance with the procedures suggested Ref. 12. Fiber content was 
determined accordance with the procedure suggested Lynn, al. (11). 

order develop visual appreciation the peats tested, their microstructure 
was examined under scanning electron microscope. examination the 
scanning electron micrographs the peat samples, given Fig. indicates 
that Middleton Peat made relatively uniformly shaped fibers with little 
amount soil grains particles between. Fond Lac Peat, the other 
hand, has higher amount soil grains and particles than the other two types 
peat. The differences fiber content are reflected the micrographs 
they were determined the laboratory. was noted that Portage Peat contains 
certain woody pieces larger than the fibers (not shown Fig. 1). 
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Available laboratory and situ techniques determine were classified 
and presented summary form Abdelhamid and Krizek (1). Within the 
framework their classification system, the technique employed this inves- 
tigation can categorized direct approach the semirigid ring type with 
continuous measurement throughout the entire consolidation history 
the peat, from its natural state with high water content consolidated 
block with relatively low water content. The test apparatus consisted 
stainless steel with strain gages mounted its wall 
monitor the time dependent history the lateral pressure. 

Tube.—Essentially, the apparatus consists honed stainless steel 
cylindrical tube (72.5-mm ID, 250 height, and thickness) that 
sealed lower end plate and piston which sealed against the cylinder 


PRESSURE SOURCE 


BELLOFRAM 


SWITCHING AND 
BALANCING UNIT STRAIN INDICATOR 


)} | | 
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wall O-ring. schematic diagram the assembled K,-test tube shown 
Fig. Vertical pressures are applied means air cylinder (bellofram) 
connected the piston. pressure difference approx kPa was found 
sufficient move the piston either down the cylinder. The piston 
friction was taken into account applying the vertical loads calibrating 
the cell. Vertical drainage was provided cast-in-place porous stone located 
the piston and the pore pressure was monitored using integrated circuit 
pressure transducer connected cast-in-place porous stone located the 
lower end plate. Starting from initial average sample thickness about 
150 mm, block peat about thick was obtained the end 
typical test. 

The side friction force between the K,-test tube and the soil specimen was 
measured determining the difference between the applied vertical load and 
the load transferred the bottom the specimen. The transferred load was 
measured load cell placed under plate below the specimen inside the 
tube which, turn, rested the loading table. The side friction force 
determinations were made for 20-30 Ottawa sand, Hydrite clay, and Middleton 
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Peat using this setup. The test results indicated that side friction increased 
with strain and pore pressure dissipation. Side friction force corresponding 
defined terms complete dissipation mea- 
sureable excess pore pressure, was linearly related the applied vertical stress 
logarithmic plot with average slope for all three soil types. The 
ratio the maximum side friction force the applied load increment per 
millimeter specimen height was, the average, 0.45-0.53 (corresponding 
0.53 and 0.66, respectively) for Ottawa sand, 0.58 for Hydrite 
clay, and 0.56 for Middleton Peat. These values fall below the values reported 
for bare steel and slightly above the values for teflon-lined rings (10); this may 
due honing and lubrication with grease the steel tube surface. 
Measurement Total Lateral hoop strains induced the 
consolidating sample against the tube wall were measured electrical resistance 
strain gages. They were 6.35-mm foil type strain gages with self-temperature 
compensation (Micromeasurements Model EA-06-250 BF-350). The gage resis- 
tance was 3.50.0 0.15% ohms and had gage factor 2.12 0.5% ohms 
24°C. They were mounted the outside surface the K,-test tube 
three levels (50 mm, mm, and 140 from the base). each level, four 
strain gages were mounted that each pair strain gages (one axial and 
one circumferential) placed 180° apart. The strain gages each level were 
connected together, thus forming full bridge arrangement with circumferential 
gages wired opposite arms the bridge. This arrangement provides bridge 
completion well temperature compensation while eliminating bending and 
torsional effects. The advantage having the four gages connected together 
that the strain sensed each strain gage added together and, thus, the 
total strain indicated the strain indicator twice the actual hoop strain 
assuming axial strain from the soil-tube friction negligible. This procedure, 
course, increases the accuracy the determined lateral pressure and, 
consequently, the coefficient. The magnitude strain was measured 
strain indicator (Baldwin-Lima Hamilton Model 120C) connection with 
switching and balancing unit (Baldwin-Lima Hamilton Model 225). The pore-water 
pressure was measured hybrid integrated-circuit type pressure transducer 
(National Semiconductors Model LX1730GB). This relatively new transducer 
offers sufficient accuracy with lower cost, smaller volume change, and greater 
flexibility over the conventional strain gage-type transducer. Furthermore, 
requires only simple and inexpensive read-out and recording systems. 
Calibration.—The strain gages have calibrated accurately order 
determine the lateral pressure exerted the soil specimen normal the wall 
the consolidation tube when vertical pressure applied. order calibrate 
the strain gages known value fluid pressure (air water) applied inside 
the tube and, since the fluid pressure equal all directions, the magnitude 
the strain measured the strain gage, corresponds the known value 
that lateral pressure. The calibration procedure indicated that all strain gages 
all levels produced linear relationships maintaining calibration factor constant 
throughout the range applied pressures (5). hysteresis effects were noticed 
during the cycling load and the strain gages all levels produced approximately 
the same calibration curve despite the difference their locations. The calibrations 
performed before and after each test showed that the electronics the 
system was stable without any significant long-term drifts. 
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The thickness the K,-test tube was chosen assure that the tube was 
sufficiently thin that the hoop deformations would accurately sensed 
the strain gages and also sufficiently thick satisfy the condition (no lateral 
deformation). Computations based the theory elasticity using the tube 
dimensions and properties, well the characteristics the strain gages, 
suggested thickness 1.31 compatible with the level sensitivity 
desired kPa lateral pressure). The tube selected had wall thickness 
about mm. The radial strain the maximum soil stress computed 
the order which significantly less than the values given 
Terzaghi (16) and Bishop (3) for the coefficient lateral earth pressure 
drop the active state. 
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order further verify that the condition satisfied, two inorganic 
soils for which the values had been reported the literature were also 
tested using the K,-test tube. Furthermore, these tests provided basis for 
the comparison the response peat soils. The two soils chosen were 
20-30 Ottawa Sand, cohesionless soil, and Hydrite PX, kaolinitic clay. 
linear relationship between applied vertical pressure and resulting lateral pressure 
obtained three different initial void ratios shown Fig. for Ottawa 
Sand. The slope these lines indicates the value for normal loading 
the corresponding void ratio. The measured values appear reasonable 
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for this sand and compare well with the values calculated using Jaky’s expression 
(Table given Ref. (8) 


This expression has been reported best represent normally loaded cohesion- 
less soils. The results the direct measurement for the kaolinitic clay 
(Hydrite PX) given Fig. 3(b). NaOH-dispersed slurry this clay (at 
initial water content 112% with concentration 0.01M NaOH solution) 
was placed the tube and allowed consolidate under increasing 
vertical stresses. The measured 0.65 for normal loading stage approxi- 
mately constant and compares well with the value 0.68 measured for similarly 
prepared kaolin using different apparatus Abdelhamid and Krizek (1) and 
with the value calculated using the Brooker and Ireland Expression (Table 
(4) 


This expression was found more appropriate for normally consolidated 
clays than Jaky’s expression. The test clay also provided check the 


reliability and stability the measurement system long-term consolidation 
tests. 


TABLE Measurements for Sand and Clay 


Determination 
from 


Effective angle 


internal Direct 
friction, measure- 


(2) 


Ottawa sand 20-30 


Kaolinitic clay hydrite 
16.9 


“Effective cohesion intercept, 


Test Procedure.—After obtaining reasonable values for sand and clay, 
peat specimens were placed the K,-test tube for measurement extruding 
them from the sampling tubes into the K,-test tube. Peat samples had the same 
diameter that the K,-test tube, since they were retrieved from the ground 
using similar tube. This made the sample preparation much easier and minimized 
the disturbance the samples. Specimens about 150-mm thickness were 
placed the tube and was expected that the final thickness after consolidation 
would about mm, corresponding total deformation about 60%. 
This final thickness would sufficient insure that the :ssure against the 
consolidation tube could reliably sensed least the iirst level gages 
which from the base the tube. All tests were performed 
constant temperature 20°C 0.5°C. general, six different quantities 
were monitored document the complete stress history the sample: time, 


K,=1- K,=0.95- 
Soil (Ref. (Ref. 
(1) (4) (5) 
0.50 0.49 0.44 
0.44 0.45 0.40 
0.41 0.43 0.38 


208 FEBRUARY 1981 GT2 


vertical deformation, vertical stress, lateral stress, pore-water pressure 
the base and effluent outflow. Loading was applied the peat specimen 
using pressure increment ratio For certain increments, the top and bottom 
drainage were initially prevented, thereby permitting the sample deform 
compressing the gas within the soil mass, thus allowing the pore-water pressure 
build its maximum value. When the highest value pore-water pressure 
was reached, and the deformation ceased, the consolidation was commenced 
opening only the top drainage line. The effluent drained from the top 
the specimen was collected burette. Each increment was left the specimen 
for period approx days which usually was sufficient complete the 
primary consolidation and considerable amount secondary compression. 
the advanced increments pressure, because the drastic decrease 
permeability, longer time was needed for the completion pore-water pressure 
dissipation and, therefore, the final pressure increments were left the specimen 
much longer period time than was the case the conventional consolidation 
tests. For example, Portage Peat, the final pressure increment (400 
kPa) was maintained the specimen for about months. The same loading 
procedure adopted for Portage Peat was followed for the other two types 
peak except, the latter two peats, loading was carried out only 400 kPa 
rather than 800 kPa. The effect overconsolidation was investigated 
unloading 50% the existing vertical pressure time and maintaining the 
specimen under that reduced pressure for period day-2 days. Before 
and after each test, the tube was calibrated evaluate the variation, 
any, strain gage response applied pressure. This procedure indicated 


that only slight changes took place, and they were within the accepted experimental 
variability. 


Test 


From the measured quantities, the total vertical and horizontal stresses, 
well the vertical strain, can determined any stage. The effluent outflow 
normalized dividing with the initial volume the specimen the 
beginning the stress increment application and referred normalized 
outflow, After excess pore-water pressure dissipates under given 
vertical stress increment, the measured vertical and horizontal stresses can 
considered effective stresses and value can computed. 

Time-dependent response during typical consolidation stage for Middleton 
Peat given Fig. this test the drainage was allowed soon the 
pressure increment was applied. Also shown this figure the ratio lateral 
vertical pressure increments, versus the logarithm time. The pressure 
ratio, gradually decreases from immediate value 0.60 with time until 
reaches value 0.32, which stabilizes. The value continues 
decrease with time even after pore-water pressure has dissipated. This behavior 
dissimilar the one observed for inorganic soils where found 
constant after complete dissipation excess pore-water pressure. Newlin (14) 
observed that lateral pressure values for compacted clay specimen stabilized 
almost exactly the same time excess pore pressure reached negligible value. 
shown Fig. the pore pressure dissipates after about 3,000 min, but 
the value does not reach its final value before 14,000 min, thus, even under 
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the same effective pressure, continues decrease. The continuing decrease 
lateral pressure indicates that considerable soil structure creep takes place. 
Apparently, the pore-water pressure associated with this process very small 
sensed the specimen boundary. The water continues flow out 
the specimen even after measurable pore pressure completely dissipated. 
This also supports the contention that pore-water pressure exists the specimen 
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FIG. Response during Typical Consolidation Stage 


even though may very small measure. The vertical strain (€) and normalized 
outflow (AV/V,) curves generally have the same shape except that the rate 
outflow slightly lower than the rate strain. For example, under this 
particular pressure increment, the deformation that can accounted for 
water outflow about mm, whereas the total deformation mm. The 
difference between the outflow and total deformation partially due the 
fact that the specimen was not saturated. However, the compressibility the 
solid phase (fibers) could also contribute this difference. The usual shape 
the primary consolidation curve observed for inorganic soils not apparent 
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compression curve. The second straight line portion which may 
termed compression (6) indefinitely continues even after the pressure 
ratio, stabilizes. some other tests, period nearly year was required 
before the measureable vertical strain ceases (6). The behavior shown Fig. 
was observed for other pressure increments with certain notable changes. 
For example, the time for complete dissipation measurable pore pressure 
increases the higher pressure increments due the drastic decrease 
permeability. The time-dependent response for the other two peat samples was 
also qualitatively similar the ones for Middleton Peat with certain 
notable differences the quantities measured. 

Void Ratio Versus Consolidation Pressure.—The curves for the three 
types peat are given Fig. Middleton and Portage Peats had relatively 
high initial void ratios, 8.0 and 7.4 respectively, whereas Fond Lac Peat, 
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FIG. 5.—Void Ratio Versus Consolidation Pressure 


which amorphous granular peat unlike the other two fibrous peats, had 
initial void ratio about 4.5. The slopes the virgin curves are almost 
identical for Middleton and Portage Peats with average compression index 
5.6. The data plotted Fig. suggests that the compression behavior 
the fibrous peat samples during virgin consolidation can divided into two 
distinct phases. The compression index decreases average value 2.5 
for these two peats after consolidation pressure 200 kPa. Fond Lac 
Peat, the other hand, exhibits constant compression index 1.5. the 
unloading stage, all peats exhibited relatively low rebound with average swell 
index about 0.3. The void ratios decreased drastically during consolidation 
and they were 2.35, 2.70, and 2.50 vertical pressure 400 kPa respectively 
for Middleton, Portage and Fond Lac Peats. 

Response.—The final values lateral pressure ratio, when stabilizes 
constant value, defined for the peat samples; stabilizes 
time beyond the defined terms the complete 
dissipation measurable excess pore pressure. Several important aspects 
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the nature and variation with stress history are manifested Fig. 
for the three peat types tested herein. The most striking observation the 
nearly linear trend the relation between the lateral and vertical effective 
stresses during the loading portions the test; this reflected the virtually 
constant values throughout the entire range loading. Furthermore, 
the values are markedly lower for peats than inorganic cohesive soils 
and the fibrous peats have lower than amorphous granular peat possibly 
due the lateral reinforcement fibers. The tensile strength fibrous peats 
known considerable (12). Adams (2) also reports relatively low values 
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for fibrous peat even though some his findings are contradiction 
with the trends reported herein, possibly due the differences the measurement 
techniques respectively used. During unloading, values become greater than 
those during loading, would expected overconsolidated soils, and the 
qualitative responses are essentially the same for the three peats with notable 
differences the quantities. However, the increase during unloading 
way comparable increase observed with inorganic soils. 
has been reported that for high values overconsolidation ratio, the value 
appears approach the value the coefficient passive earth 
pressure and generally has value greater than one (4,14). For the three peats 
tested, the value always less than one overconsolidation ratio 
shown Fig. Despite the possibility tube wall contraction resulting 
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higher values during unloading than really exist the soil, the qualitative 
trends are probably quite reasonable. Schmidt (15) has proposed the following 
formula relate (during rebound) and (during normal loading): 


which OCR the overconsolidation ratio; and coefficient. The average 
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Middleton Peat 


values were found 0.18 for Middleton and Fond Lac Peats, 
0.09 for Portage Peat, and 0.19 for Hydrite clay. This equation represents 
the measured function OCR reasonably well. 


Comparison 


The foregoing direct measurements are compared utilizing various 
indirect approaches conjunction with the strength parameters from shear 
testing program the same samples used determine values K,. series 
anisotropically consolidated undrained triaxial compression tests were con- 
ducted different consolidation pressures determine the shear strength 
characteristics the peat samples. The stress-strain response and the induced 
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pore pressures are shown Fig. along with the total and effective failure 
envelopes for sample Middleton Peat. The general patterns indicated 
Fig. were observed for the other peats. The induced pore pressures during 
undrained shear were found relatively low and the pore pressure parameter, 
A,, was generally less than 0.5 for all peats. The shear strength mainly 
frictional with exceptionally high values the effective angle internal friction, 
(50°-57°) compared those inorganic soils. Similarly high values 
were also reported Adams (2) indicating that the shear strength peat 
greatly increases consolidation. 

number expressions, mostly based semiempirical considerations, have 
been suggested evaluate from measured values Value obtained 
from each these expressions, which were developed basically for inorganic 
soils, for the measured values are given Table along with the measured 
values K,. the calculated values with those measured directly 


TABLE Determination from Shear Strength Parameters 


Most Middleton Fon Lac 
appropriate Peat Portage Peat Peat 
(1) (2) (3) (4) (5) 
1-1.2 sin (Ref. 17) Cohesive 


8 8 


(Ref. 7) 
Direct Measurement 


Cohesionless 


suggests that none these equations are satisfactory for peat soils. Such 
discrepancy could attributed many factors including the argument that 
the empirical equations are valid only for certain range values and 
for certain soil types. Lambe (9) suggests that the original theoretical derivation 
given Jaky the form 


valid for values from 20°-45°. apparent that not adequate, 
single parameter, reflect the lateral load transfer mechanisms due 
the different nature the microstructure peat soils, particular, due 
the presence fibers. Despite the general inability express the behavior 
the peats function certain fundamental parameters, i.e., friction angle, 
fiber content, etc., possible distinguish broadly between the behavior 
fibrous and amorphous granular peats. shou!d pointed out that Middleton 
and Portage Peats both are classified fibrous, whereas Fond Lac Peat 
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amorphous granular, according the classification MacFarlane and 
Radforth (13). is, therefore, suggested that low value about 0.30 
associated with fibrous peats, whereas amorphous granular peats have higher 
about 0.53. This sort broad distinction between fibrous and amorphous 
granular peat groups was also observed investigation compression behavior 
involving large number tests and four different peats (6) and supports the 
contention that this classification meaningful for engineering purposes. 


view the experimental results and the analyses presented this 
investigation, the following conclusions can advanced: 


The tube instrumented measure lateral stresses and pore water 
pressures during one-dimensional consolidation provides unique experimental 
capability for investigating the behavior wide variety soils from 
sands clays peats. 

The ratio the horizontal stress vertical stress the time stabilizes 
constant value (usually some time after the consolidated equilibrium 
reached) defined the coefficient lateral earth pressure at-rest, for 
peat soils; constant throughout the range virgin loading. 

Significant difference observed between the values peats and 
clays, with the former being much lower due the presence fibers; furthermore, 
peat type exerts significant influence the measured value during loading, 
which 0.53 for amorphous granular peat and about 0.30, markedly low 
value, for fibrous peats. 

The relationship between the horizontal and vertical effective stresses during 
unloading nonlinear, and the value is, therefore, not constant; however, 
the increases during unloading not significant those observed 
for clays. 

The effective friction angle peats much higher than that for inorganic 
soils (about 50° for amorphous granular and 53°-57° for fibrous peats); the 
use the effective friction angle determine indirectly from empirical 
expressions was found give unacceptable values when compared with 
those measured the laboratory. The use the classical soil mechanics strength 
theory explain the behavior peats limited due the presence fibers. 

Classification peats broadly into the categories fibrous and amorphous 
granular peats appears meaningful for engineering purposes; this supported 
the behavior reported herein well the compression behavior 
reported previously (6). 


The value soil determines the magnitude situ shear stresses 
soil deposit. The relatively low values for peat soils implies the 
presence higher situ shear stresses for given effective vertical stress 
than those clay deposits. important recognize this both with respect 
the stability problems the field involving peats well the usefulness 
isotropically consolidated triaxial compression testing peats the laboratory. 
The effect fibers observed for fibrous peat and certain industrial 
sludges (such papermill sludge) leads one conclude that normal soil mechanics 


216 FEBRUARY 1981 GT2 


parameters cannot used for such materials. The writers further recognize 
that environmental changes could have effect fiber quality and strength 
and consequently the parameters examined could affected over time with 
certain environmental changes. 
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The following symbols are used this paper: 


pore pressure coefficient; 

compression index; 

swell index; 

effective cohesion intercept; 

void ratio; 

coefficient; 

ratio lateral vertical pressure; 
coefficient earth pressure 
coefficient earth pressure at-rest during normal loading; 
coefficient earth pressure at-rest during rebound; 
overconsolidation ratio; 

excess pore pressure; 

normalized outflow; 

axial strain; 

major principal stress; 

vertical consolidation pressure; 

minor principal stress; 

effective horizontal stress; 

effective vertical stress; 

consolidated undrained angle friction; and 
effective angle friction. 


e = 
K = 
K onc) = 
OCR 
AV/V, 
€ => 
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STATISTICAL ANALYSIS MARINE CLAY 


Hin Fatt ASCE 
and 


INTRODUCTION 


Extensive deposits marine clay are found around the regions 
Singapore where several projects are currently progress. The wide variations 
observed physical properties such marine clays make statistical approach 
useful tool. 

Skempton (4) suggested relation c/o 0.11 0.0037 (w, w,) which 
effective overburden pressure; liquid limit, and plastic limit. 
Osterman (3) also proposed similar linear relationship between the plasticity 
index and c/o for marine clays. Koutsoftas and Fischer (2) reported graphically 
the relation among the undrained shear strength, overconsolidation ratio, and 
the effective consolidation stress for marine clays along the New Jersey coast 
the United States. Balasubramaniam, al., (1) expressed the undrained strength 
Bangkok clays terms natural water content, unit weight, and effective 
overburden pressure. 

the present investigation, the main objective find dimensionless 
regression equation for the normalized undrained strength Singapore marine 
clays terms other normalized properties easily determined the laboratory. 

Fig. summarizes the sequence analytical operations and the reasons for 
the same. 


Preliminary analyses the data obtained from bore logs indicate that the 
mean values the same soil property for marine clay deposits the various 
sites are significantly different. The soil parameters can studied normalizing 
each soil property according the equation 
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which ith observation the soil property, standard deviation 
Y,. Normalization the data from the various sites will give rise normalized 
data sets from each site. 

Since the normalized variable has expected zero mean and unit variance 
irrespective the site, conceivable that the sets for each site are 
statistically similar. other words, possible that for each site can 
described the same probability law. 

The assumption first made that the sets describing any soil property 
for the various sites are statistically similar. The sets from the different 


FLOW CHART 


H-test on raw data 


To ascertain whether raw data from various sites 
belong to the same population 


Normalisation of raw data 
It is possible that normalised variates with zero mean 
and unit variance are from the same general population 


Combined normalised data 


If the above postulate is true, combining the normalised data 
is valid and this results in a larger sample for statistical 
analysis 


Stepwise regression analysis 


To check if postulate is acceptable, a multiple regression equation 
Obtained through stepwise analysis is made. The residuals are then 
determined. Kruskal-Wallis H-test is required to determine whether 
the residuals are from the same population. If H-test is affirmative, 
then normalisation and grouping of normalised data is acceptable. 


Randomness test on residuals 


A basic requirement in the H-test is that the 
residuals are random variables 


H-test on residuals 


To determine whether the residuals from various 
sites belong to the same population 


Frequency distribution 
Necessary for sequential analysis 


Sequential analysis 


To check validity of regression 
analysis at a new site 


confidence intervals 
FIG. 1.—Flow Chart 
sites are then grouped together form larger set which 


The set residuals the normalized cohesion, c*, for given site, 
obtained according the following equation 


which predicted value obtained from the stepwise regression 
equation. The set residuals are analyzed ascertain: (1) they are random; 
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(2) their frequency distribution; and (3) they are independent the normalized 
variates the regression equation. 

The Wald-Wolfowitz Runs test being used ascertain the randomness 
the set residuals given site with the residuals taken order 
occurrence with increasing effective stress. The test randomness has 
applied every set residuals from all the sites before the application 
the Kruskal-Wallis H-test for deciding whether the residuals the normalized 
data sets from the different sites come from identical populations. The amount 
data collected new site where marine clay concerned generally 
quite limited, especially the area under development not extensive. these 
data values are normalized and can grouped with normalized values the 
sites from which the regression equation obtained, then the dimensionless 
stepwise regression equation derived earlier may also apply the new site. 
The sequential probability ratio (SPR) test may used ascertain the expected 
number samples required establish whether grouping justified not. 

Let the set residuals for the new site, and P,(r) and P,(r) 
the probability density functions with unit variance and means and 
respectively. desired test the hypothesis that the observed set 
came from the distribution P,(r) against the possibility that came from 
another distribution P,(r). The proposed sequential test employs the likelihood 
ratios sequentially. presupposes that constitutes random sample; i.e., 
the r,s are independent and identically distributed. The procedure entails the 
definition the likelihood ratio 


i=! 
and the sequential computation computing these ratios. each step the 
computation, necessary determine the ratio lies within the interval 
for preassigned Type and Type errors and respectively). 
The value the probability rejecting when true, and the 
probability accepting when false. The choice and fixes 
since 


The test continues beyond the jth step long and stops 
the mth step when A,, k,. The hypothesis rejected 

For the case which P,(r) and P,(r) are both Gaussian distributed with 
means and and standard deviation 


formean 
i=l 


and formean 
i=l 
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Resutts 


Data from three different geographical sites, namely Changi, Buona Vista, 
and East Lagoon Singapore are analyzed. Since different testing methods 
the soil samples yield different results, only those data covered the undrained 


TABLE 1.—Results Stepwise Regression Analysis Based Linear Combinations 
Independent Variables 


Standard 
error 
estimate 


Computed 
statistic 


(3) 


Variable 


(1) 


(a) Step 


44.7 


R(w*) 27.8 3.92 
3.92 
20.9 


4.6 


Step 


3.65 


0.61 3.00 
2.05 3.00 
20.57 3.00 
39.07 3.00 


25.89 
40.45 


(c) Step 
0.32 
0.28 
0.43 


number independent variables. 

Note: Regression equation: 0.480* 0.2w*; multiple regression coefficient 
0.62. 


triaxial test are considered. The normalized engineering properties considered 
are effective stress, cohesion, c*; bulk density, y*; water content, w*; 
liquid limit, plastic limit, w}; liquid index, and plasticity index, 
Bore log data marine clay without information all the aforementioned 
properties are not considered. Ambiguous data that accrued from tests performed 
samples near boundaries stratification are also discarded. summary, 
there are 77, 38, and data sets from Changi, Buona Vista, and East Lagoon, 


with 
| 
| 
(4) 
0.80 3.92 
0.79 3.00 
0.52 2.60 
0.46 2.60 
0.69 2.60 
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respectively. The maximum depths which samples have been taken about 
while 80% the samples are from depths less than below ground 
level. Overconsolidated clays have not been included. For all sites, the bulk 
density exhibits the lowest variability about the mean value. Consequently, the 
use the average value the estimation the effective stress any prescribed 
depth due the overlying surcharge will not incur serious errors. 

The results the stepwise regression analysis are summarized Table 
and the analysis based equations expressing c*, the dependent variable, 
terms linear combinations the independent normalized variables. 


c 


WHERE 


LEGEND 
95%. CONFIDENCE INTERVAL 
FOR C 
CONFIDENCE INTERVAL 
FOR TRUE MEAN 
CHANG! 
BOUNA VISTA 
EAST LAGOON 
NEW SITE 


FIG. Regression Equation 
The stepwise regression equations (SRE) are found 


with multiple regression coefficients 0.62 and 0.61, respectively. Eq. 
preferred, does not involve the determination Atterberg limits. The 
graph the regression equation with the confidence intervals are shown 
Fig. with the normalized data from the various sites plotted for comparison. 
For each site, the residuals the normalized variates about the SRE are 
obtained and the Kruskal-Wallis H-test performed for the three sets residuals 


| 
a 
2 
40 
30 
4 
. 
sa 6 
o* 
-30 Ne 10 20 30 x 
a 
oN x 
‘\ 
-3.0 
\ 
-40 
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from Changi, Buona Vista, and East Lagoon. The computed H-statistic 0.62, 
which less than the critical value 5.99 the level significance. 
This test result indicates that the grouping the three sets residuals 
justified. the Two-Way Contingency Table tests, the residuals are ascertained 


TABLE 2.—Summary Sequential Probability Ratio Test Marine Clay Parameters 
New Site 


per square meter per square meter 
(2)° (3)° 


w, asa 
percentage 


“Mean 14.37; standard deviation 2.49. 
111.47; standard deviation 49.47. 
“Mean 61.50; standard deviation 4.32. 


TABLE 3.—Summary Sequential Probability Ratio Test Marine Clay Parameters 
New Site 


3.277 
2.557 

4.293 

2.979 

2.277 
20.28 (Stop) 


TABLE 4.—Summary Sequential Probability Ratio Test Marine Clay Parameters 
New Site 


0.829 
2.887 
4.675 
18.305 

65.11 (Stop) 


illustration the application the SRE for estimating the cohesion 
values marine clay new locality, the residuals from the SRE are determined 
and the SPR test used determine the normalized values the new 


(1) 
10.5 42.38 
17.2 150.09 
12.9 169.71 
(1) (4) (5) 
0.584 0.748 0.780 
1.678 
0.675 0.865 0.694 
0.600 0.769 0.764 
Data 
(1) (4) (5) 
0.584 0.748 3.483 
1.619 
0.675 0.865 3.915 
0.600 0.769 3.557 
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site can incorporated into the total grouping. The results such analysis 
based Gaussian assumption for residuals given Tables Table 
the set the first six data taken. Tables and the residuals are 
normalized using 0.78 and the likelihood ratios, \,, are computed. 

These results are obtained the basis limiting the Type and Type 
seen that least five six samples are required ascertain the normalized 
variates marine clay the chosen new site can grouped into the SRE; 
five samples are required for the case which postulated have 
and variance, and six for the case which postulated 


have mean and variance Upon substitution the means and variances 


This approach especially use there are limited data for the new site where 
lack adequate data can render regression analysis meaningless. 


Marine clay properties from the various geographical areas Singapore could 
not grouped together for statistical analysis their raw form. This necessitated 
normalization the areal level and the application the Kruskal-Wallis H-test 
the residuals showed that the normalized variates can grouped for further 
statistical analyses. The final forms the stepwise regression analysis which 
involved the independent variables o*, y*, are 


which the multiple correlation coefficients are 0.62 and 0.61, respectively. 
The validity these equations has not been tested for overconsolidated and 
cemented clays. 

For the Type and Type errors, the SPR test shows that small 
number samples sufficient infer whether the normalized variates 
the new site can grouped with the total sample used earlier for regression 
analysis. This test, affirmative, also permits the use the stepwise regression 


equation for reduction simple predictive equation for cohesion the new 
site. 
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PLASTICITY APPROACH SLOPE 
STABILITY PROBLEMS 


Ryszard 


INTRODUCTION 


slope stability analysis, various limit equilibrium methods are applied, e.g., 
method (4), Bishop’s method (1), Janbu’s method (6), Kenney’s method 
(8), the Morgenstern-Price method (11), Madej’s method (10), and others. All 
these methods are generally called The safety factor 
obtained from the conditions global equilibrium forces acting the 
vertical rigid slices (see Fig. 1). Thus, the methods are similar the limit 
analysis method, which based the static equilibrium the field forces 
being associated with the failure mechanism (see Ref. 5). However, application 
the slices methods requires various additional assumptions concerning the 
force distribution interfaces slices, since the stability problems are statically 
indeterminate. 

The application the upper bound theorem the determination the safety 
factor examined this paper. The upper bound the safety factor compared 
with those obtained conventional limit equilibrium methods. The accuracy 
limit equilibrium methods also estimated. 


classical measure the margin slope safety the safety factor, 
defined relation between shear strength and mobilized shear stress; i.e. 
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which shear strength; mobilized shear stress; normal stress; 
tan friction; cohesion; tan mobilized friction; and mobilized 
cohesion. The factor assumed have the same values along both shear 
surfaces and interfaces slices. 

Assuming the associated flow rule, obvious inequalities follow from limit 
theorems (see Ref. 2); thus 


which the real value the safety factor; and and the lower 
and upper bounds, respectively. 

The upper bound the safety factor, F,, obtained assuming that along 
shear surfaces and along interfaces slices there occurs the limit state 
stress given If, the slices methods (1,4,6,8,10,11), consent 
the existence along interfaces slices the stress state below the condition 


FIG. 1.—System Forces 


then the determined value the safety factor, lies below the 
upper bound: 


This inequality expresses the theorem the reduced upper bound. 

Eq. may reasoned the following way. Let assume that along interfaces 
slices there thin layer material obeying the yield condition 
which lies totally inside the limit surface corresponding the condition 


The upper bound theorem valid but since the 
result Eq. 


When applying the limit analysis method the upper bound can obtained 
the following methods: (1) comparing the sum the total work external 
forces and total work body forces with the total dissipation energy (see 
Refs. 7); (2) making use the equilibrium conditions the field 
force associated with the assumed kinematically-admissible failure mechanism 
(see Refs. 9). 


T 
a 
bd 
| Toi Ni 
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Both approaches are equivalent. can observed that the equation the 
total work balance may treated equation that expresses the condition 
global equilibrium the system consisting rigid blocks. Also, the local 
equilibrium each block may performed, provided that aiong the division 
lines the shear stresses are limit value. 

Let consider now the system forces associated with the failure mechanism 
which the sliding mass divided into vertical rigid slices (see Fig. 1). When 
introduce the symbols: AE, E,,, E,; and AX, X,,, the 
equilibrium equations may written 


and express the horizontal and vertical equilibrium; Eq. expresses 
the moments equilibrium about the center point the slice base. 


FIG. 2.—Failure Mechanism 


Assuming the forces along the interfaces slices and base the same 
those for the limit state; i.e. 


of F, ’ F, ’ F, (5) 


and making use Eqs. and 4b, may write the condition internal 
force equilibrium for the whole sliding mass, the following form 


which 
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k 


k 
From these equations may directly determine the value the upper bound 
the safety factor. The calculation the increments AX, unnecessary, 


FIG. 3.—Force Polygon 


F, 
From the aforementioned, follows that when the condition 
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For the slope high (see Fig. 2), with 1:2 inclination and soil parameters 
2,17. The piecewise linear and circumscribed the circular arc slip surface, 
were assumed the shear surface bounding the sliding mass. 

The satisfaction equilibrium equations geometrical interpretation corre- 
sponds the building the closed polygon forces which the case considered 
may obtained building successively the polygons for slices and 
(see Fig. 3). should added that the assumption limit values 


interactions along interfaces slices the stability problem statically 
determinate. 


FIG. Hodograf 


When constructing the hodograf velocity (see Fig. and making use 
the equation total work balance, again arrive the value 2,17. 
The assumed failure mechanism consists shearing along the shear surface 
the sliding mass and along the vertical surfaces dividing the sliding 
mass into slices (see the translation scheme Fig. 2). Internal dissipation was 
calculated from the formula 


and total work body forces was given 


0 
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which velocity slice; and determines the velocity discontinuity 
along the appropriate vertical line separating the neighboring slices. 


The calculation example was taken from Ref. (also see Ref. 11), which 
the safety factors are given determined some slices methods. The application 
the correction Kenney’s method (10) resulted 2,045; the correction 
Janbu’s method gave 2.050, whereas Madej’s method, 2.048. 

All the methods allow the construction the associated force field 
equilibrium; this field, however, the state forces interfaces slices 
lies below the limit state. Thus, obtain the bounds safety factor according 

Assuming, additionally, that there are forces along the interfaces slices 
(as Fellenius’ method), obtain 1.95 from global equilibrium conditions. 
this case, local conditions slices equilibrium are satisfied for each slice 
separately. So, may state that all the accurate methods slices, provided 
that the yield condition not violated, will give factor bounds the interval 


The upper bound approach limit analysis presented the paper allows 
the upper bound rapidly obtained means analytical graphical 
method. The use other methods mentioned for calculations most cases 
troublesome because requires the application electronic calculation 


technique, and the increase calculation accurracy for any the methods 
not proportional the increased expenditure work and energy. 
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manuscripts. 


DISCUSSION 


SHEAR AND LIQUEFACTION 
Discussion Yoshiaki Yoshimi,’ ASCE and Kohji 


The authors have made valuable contribution resolving some the 
conflicting views the influence initial static shear stress the behavior 
sand subjected cyclic shear under constant volume conditions. 

The writers admit that the simple shear test apparatus used the authors 
quite suitable for carrying out the constant volume cyclic loading tests. The 


ability the apparatus apply large shear strains large shear stress ratios 
may attributed the fact that displacements rather than stresses are prescribed 
along the boundary surfaces the specimen. other words, the specimen 
forced deform into parallelepiped with well-defined average shear 
strain, the expense nonuniform stresses along the boundaries. 

the other hand, the ring torsion shear apparatus developed the first 
writer and Oh-Oka (8,10) was designed reduce the nonuniformity the 
boundary shear stresses making the specimen endless the direction 
shear and making the height the specimen proportional the radius 
shown Fig. 11. the authors point out, the ring torsion apparatus does 
impose nonuniform stresses and strains during consolidation. The first writer 
and Oh-Oka have shown, however, that such nonuniformity was considered 
insignificant when the geometry the sample was shown Fig. (8). 

order facilitate comparison, the test results the authors and those 
the first writer and Oh-Oka (8) have been replotted Fig. terms 
normalized quantities. The ordinate shows the cyclic shear stress ratio, 
required cause 2%, 10% shear strain stress cycles any value 
static cyclic shear stress ratio, normalized with respect that 
under condition total shear stress reversal, indicated above 
total, partial, and reversal shear stress, respectively. The fact that 
the curves for 37% lie below those for 50% shows that the test 
data the first writer and Oh-Oka are compatible with the authors’. The 
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Tokyo Inst. Tech., 4259 Nagatsuta, Midori-ku, Yokohama 227, Japan. 

Prof., Tokyo Inst. Tech., 4259 Nagatsuta, Midori-ku, Yokohama 227, Japan. 


GT2 237 


238 FEBRUARY 1981 GT2 


dimensionless relationship Fig. confirms the authors’ important conclusion 
that the cyclic shear strength sand may 


increase, decrease remain unaltered, depending the relative density 


the sample, the shear strain level interest and the magnitude 
initial static shear stress. 


Although Fig. represents constant amplitude cyclic shear tests which 
was held constant during each test, the value varies during 


FIG. 11.—Shape Specimen for Ring Torsion Shear Test First Writer and Oh-Oka 
(8) 


Total shear stress reversal 
Partial stress reversal 
shear stress 


D-= 


Usual range seismic stresses 


FIG. Test Results Authors and First Writer and Oh-Oka 
(8) 


earthquake because varies while remains constant. Therefore, the 
concept equivalent uniform cyclic shear stress (9) that has been applied 
the completely reversed shear stresses might require certain modification 
when static shear stress present. Nevertheless, one takes 65% the maximum 
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value irregular dynamic shear stress the amplitude equivalent uniform 
cyclic shear stress, the ratio embankments and level ground 
near structures during strong earthquakes seldom exceeds unity, i.e., shear stress 
reversal usually occurs indicated the bottom Fig. 12. 
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